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NOMENCLATURE 


Es  /  '  T 

A  Nondimensional  parameter,  -3-  [1  -  — ^ 

RT, 


A’  Nondimensional  parameter,  [-3 - 1/2|  [l  - 


RT. 


A^  Are?  of  burning  surface 
A£  Throat  area  of  sonic  nozzle 
a  Acoustic  wave  velocity 

Nondimensional  parameter  (see  Eq.  4.9) 

Bhfc  Nondimensional  parameter  (see  Eq.  4.10) 

C  ,  C  Specific  heat  of  solid  and  gas,  respectively 

b  p 

c*  Characteristic  velocity  of  exhaust  gases 
D  Diameter 

E,.  Ej  Activation  energy  of  the  solid  and  the  flame,  respectively 


*  _ 

j  Temperature  sensitivity  of  propellant,  r— 

r  oT 
T  o 

K  A  constant,  <_  1 

L  Length  of  T-burner 

L*  Characteristic  length  of  chamber  (free  volume  divided  by  nozzle 
throat  area) 

l g  Distance  between  entropy  extremes 

m  Mass  burning  rate 


v 


NWC  TP  4478 


n  Pressure  exponent  in  steady  state  burning  rate  law,  r  “  bp 

p  Pressure 

"p  Mean  or  time-average  pressure 

p'  Pressure  perturbation 

p^  Reference  plane  pressure 

R  Gas  constant 

Re()  Real  part  of  a  complex  expression 

r'  Burning  rate  perturbation 

s  Entropy 

sq  Reference  entropy 

T  Temperature 

T  Mean  or  time  average  temperature 

T'  Temperature  perturbation 

Tj  Adiabstic  flame  temperature 

Tq  Initial  temperature  of  the  propellant 

Reference  plane  temperature 
Tg  Surface  temperature  of  burning  propellant 

T  Period  of  acoustic  oscillation 

V  Volume 

V  Velocity  amplitude 

v  Velocity  perturbation 

v  Mean  or  time-average  velocity 

v^  Mean  gas  velocity 

X  Displacement 

V  Acoustic  admittance 


vi 


NVIC  TP  4478 


a  Exponential  growth  constant 

at  Thermal  diffusivity  of  aolid 

s  Ratio  of  maximum  presaure  perturbation  to  mean  pressure 

X  Characteristic  root  of  the  aolid  phase  energy  equation  (see 
Eq.  4.7  and  wavelength  in  Section  5) 

u/e  Ratio  of  masa  perturbation  to  pressure  perturbation.  Also 
designated  the  propellant  response  function. 

p  Density 

o/e  Ratio  of  density  perturbation  to  preaaure  perturbation 
t  Characteriatlc  time  to  burn  through  a  particle 
T(.jJ  Residence  time  of  the  chamber 

“t 

t  Thermal  wave  relaxation  time,  — 


♦  Phase  of  combuation  relative  to  pressure 

i l>  Phase  between  preaaure  and  emergence  of  element  of  gaa  from  the 

combuation  zone 

“t  L 

fl  Non-dimensional  frequency,  ~zr~ 

] r L 

u  Frequency  in  radiana/aecond 


SUPERSCRIPTS 


—  Bar  over  a  aymbol  (e.g.,  r)  refera  to  a  time  averaged  value 


SUBSCRIPTS 


b  Burning  surface 
g  Gas 


o  Total 


c  Chamber 
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INTRODUCTION 


1.1.  CONCEPT  OF  STUDIES 

The  research  work  described  in  this  report  Is  part  of  a  continuing 
program  (Ref.  1)  of  inveatigation  of  the  processes  Involved  In  low  fre¬ 
quency  acoustic  and  nonacoustic  type  oscillatory  behavior  In  solid  pro¬ 
pellant  rocket  motors  and  In  various  aspects  of  steady-state  combustion 
and  Ingredient  decomposition  believed  to  be  related  to  solid  propellant 
combustion  Instability.  The  present  program  was  Initiated  at  a  time  when 
there  were  two  concurrent  trends  In  the  development  of  solid  propellant 
rocket  motors.  The  first  of  these  was  the  development  of  very  largo 
solid  propellant  booster  rockets,  whose  dimensions  were  such  as  to  permit 
oscillatory  behavior  In  acoustic  modes  at  frequencies  down  to  about  10  cps. 
The  second  was  the  growing  Interest  In  motors  Intended  for  space  appli¬ 
cations  which  possessed  very  high  volumetric  loading  density  and  were  de¬ 
signed  to  operate  at  relatively  low  pressures;  both  of  these  features  are 
Important  for  favorable  mass  ratios  required  for  upper  stage  applications. 
Some  motors  of  this  type  were  observed  to  exhibit  oscillatory  behavior  in 
a  low  frequency,  nonacoustic  mode,  particularly  In  the  frequency  range  of 
10  to  200  cps  (Ref.  2-4).  Considering  the  extremely  high  cost  of  Individ¬ 
ual  firings  of  the  very  large  rocket  motors  under  development  and  the 
cost  of  trial  and  error  corrections  should  low  frequency  combustion  Insta¬ 
bility  occur.  It  was  advantageous  to  develop  laboratory  scale  tests  which 
would  permit  qualitative  evaluation  of  the  susceptibility  of  propellants 
to  low  frequency  combustion  Instability  and  to  determine  those  attributes 
of  the  combustion  responsible  for  Instability. 

A  review  of  the  general  character  of  solid  propellant  combustion  and 
of  oscillatory  combustion  were  presented  in  the  previous  report  (Ref.  1); 
that  review  should  be  considered  as  background  Information  for  the  work 
reported  here. 


1.2.  PUN  OF  THE  REPORT 

This  report  summarizes  work  that  has  been  done  on  this  program  from 
April  through  September,  1967.  Each  section  deals  with  a  specific  aspect 
of  the  program  and  in  some  cases  is  complete  in  Itself.  Although  the 
sections  arc  independent,  they  are  mutually  complementary;  results  in  one 
area  often  suggest  Investigations  that  should  he  pursued  in  others. 

Section  2  describes  the  continued  investigations  of  decomposition 
and  deflagration  of  large  single  ammonium  perchlorate  (AP)  crystals.  The 
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work  Included  Investigations  with  MnO^doped  crystals  and  combustion  of 
crystals  In  low  pressure  methane. 

Section  3  reports  the  Improvements  In  the  apparatus  and  techniques 
for  photographing  propellant  samples  In  a  window  bomb  and  also  the  prog¬ 
ress  In  understanding  the  detailed  aspects  of  aluminum  accumulation  in 
"real  propellants". 

Section  4  covers  the  work,  both  experimental  and  theoretical,  on  non- 
acoustic  combustion  instability.  Experimental  data  was  reasonably  well 
correlated  by  the  models  of  Denison  and  Baum  and  of  Hart,  Farrell,  and 
Cantrell.  A  phenomenological  concept  la  presented  which  partially  ex¬ 
plains  the  deviation  of  experimental  data  for  several  blmodal  oxidizer 
propellants  from  the  one-dlmenslonal  thermal  theory. 

Section  5  deals  with  the  work  on  low  frequency  acoustic  Instability 
with  particular  emphasis  on  the  dynamics  of  the  combustion  zone  as  In¬ 
ferred  from  streak  photographs  and  gas  temperature  measurements  taken  dur¬ 
ing  oscillatory  combustion. 

Section  6  is  a  brief  summary  of  the  more  significant  observations 
made  during  this  reporting  period  and  their  relation  to  previous  work  on 
the  program. 

Appendix  A  consists  of  burning  rate  data  and  tsbalar  experimental 
data  obtained  on  nonacoustic  combustion  Instability.  Appendix  B  presents 
the  results  of  testing  some  of  the  NASA  big  booster  candidate  propellants 
in  both  the  acoustic  and  nonacoustic  laboratory  burners  snd  discusses 
their  significance. 
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2.  AMMONIUM  PERCHLORATE  DECOMPOSITION 
AND  OEFLAGRATION 


2.1.  INTRODUCTION 

Ammonium  perchlorate  (AP)  Is  the  oxidizer  and  principal  Ingredient 
In  most  modern  composite  propellants  and  Its  behavior  in  the  combustion 
wave  la  an  essential  aspect  of  the  combustion  as  a  whole.  The  knowledge 
of  the  major  aspects  of  the  low  temperature  decompoaition  and  high  tem¬ 
perature  deflagration  of  AP  was  summarized  In  the  previous  report 
{Ref.  1).  Work  has  continued  in  both  these  a^eas  of  AP  behavior.  In  the 
thermal  decomposition  ares,  single  crystals  of  pure  AP  were  partially  de¬ 
composed  and  subsequently  studied  with  a  scanning  electron  microscope  and 
the  interfacial  decomposition  of  MnO^-doped  AP  crystala  was  investigated. 
In  the  area  of  AP  deflagration,  further  studies  were  conducted  to  eluci¬ 
date  the  extent  of  the  occurrence  of  a  melt  on  the  surface  of  pure  AP 
crystala  undergoing  self-deflagration.  Also,  deflagration  studies  have 
been  started  on  doped  AP  single  crystals. 


2.2.  LOW  TEMPERATURE  DECOMPOSITION  EXPERIMENTS1 

2.2.1.  Surface  Structure  of  Partially  Oecompoaed  AP  Crystals 

To  Rain  additional  insight  into  the  low  temperature  decomposition  of 
AP,  s  series  of  experiments  were  undertaken  to  systematically  study  the 
surface  features  of  partially  decomposed  AP  crystals.2  The  isothermal 
decomposition  of  single  crystala  of  AP  was  monitored  with  a  hot  ataRe 
microscope  and  stopped  at  a  predetermined  stage.  The  samples  were  then 
vacuum  coated  with  an  Au/Pd  alloy  and  kept  under  vacuum  until  they  could 
be  photographed  with  the  scanning  electron  microscope  (SEM).  Several 
serlea  of  SEM  photographs  are  shown  of  AP  samples  decomposed  under  various 
conditions.  Each  figure  contalna  pictures  of  four  different  magnifica¬ 
tions  to  allow  the  study  of  details  as  well  aa  bigger  areas  of  a  crystal 
surface.  Rectangular-shaped  holea  are  characteristic  for  decomposition 


1 

This  work  waa  supported  by  Naval  Ordnance  Systems  Command  ORO-033  129/ 

200  1/R001-06-01,  Problem  Assignment  15  and  is  reported  here  for  continuity. 

2 

Both  the  pure  and  doped  AP  crystals  uaed  in  the  studies  In  this  section 
were  grown  by  Or.  W.  R.  McBride  of  the  Research  Department,  Michelson 
Laboratory,  NWC,  China  Lake,  California.  A  description  of  the  techniques 
is  given  in  Ref.  1. 
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sices  on  m-faces  whereas  rhombohedral-shaped  holes  sre  characteristic 
for  the  decomposition  sites  on  c-faces.  Figure  2.1  shows  the  m-fsce  of 
a  crystal  which  was  decomposed  Isothsrmally  at  226°C  for  40  minutes  at 
an  sir  pressure  of  20  mm  Hg.  The  characteristic  feature  Is  the  rectan- 
gulsr  structure  of  the  holes  which  make  up  the  decomposition  nucleus. 

The  holes  have  the  shape  of  the  m-face  with  edges  running  parallel  to 
the  edges  of  the  m-face.  They  are  arranged  In  a  disk-like  pattern  and 
eventually  develop  Into  s  porous  residue.  Furthermore,  the  decomposition 
seems  to  occur  In  lsyers  with  the  size  of  the  holes  becoming  smaller  and 
smaller  with  depth.  There  is  a  second  type  of  decomposition  site  which 
generally  does  not  develop  into  growing  nuclei  composed  of  holes.  They 
sppenr  Immediately  after  reaching  the  decomposition  temperature  snd  were 
not  observed  under  atmospheric  pressure. 


— i  -H  i—  ~r***T~  h  •*'*  i — 


—| -el—  H  •<*  I—  — i*rj—  —)  •<*  I— 

FIG.  2.1.  Decomposition  Site  on  FIG.  2.2.  Decomposition  Site 

m-Face  of  Pure  AP  After  40  Minutes  on  C-Face  of  Pure  AP  After  60 

Decomposition  at  226°C  and  20mm  Hg.  Minutes  Decomposition  at  226°C 

and  20mm  Hg. 


Figure  2.2  shows  part  of  a  C-face  on  a  crystal  which  was  decomposed 
for  60  minutes  at  226°C  and  20  mm  Hg  of  air.  Analogous  to  m-faces  there 
sre  two  types  of  decomposition  sites  which  have  the  shape  of  rhombohe- 
drons.  One  type  consists  of  holes  arranged  In  a  streak-like  pattern. 

The  streaks  grow  faster  along  the  b-axls  than  along  the  s-sxls  and,  after 
overlapping,  constitute  the  porous  residue  layer.  The  other  type  com¬ 
prises  very  flat  rhombohedral  holes  which  are  spread  all  over  the  C-faces 
and  appear  only  under  low  pressure.  Both  the  streaks  and  the  long  diag¬ 
onals  of  both  typea  of  decomposition  sites  are  parallel  to  the  b-axls  of 
the  crystal.  If  the  decomposition  is  carried  on  for  a  longer  time  the 
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surface  becomes  less  regular,  although  some  of  the  rectsngular  structure 
Is  still  preserved,  as  can  be  seen  in  Fig.  2.3,  which  shows  part  of  an 
tn-fsce  after  140  minutes  decomposition  at  226°C  and  a  pressure  of  20  mm 
Hg  sir.  There  is  s  very  different  type  of  structure  at  the  very  top  of 
the  decomposition  site  in  Fig.  2.3.  It  appears  to  consist  of  clusters 
of  very  small  AP  particles.  The  clusters  may  be  products  of  a  complex 
dissociation-desorption-adsorption  process  or  may  be  residues  of  advanced 
sublimation  of  the  porous  lsyer. 


FIG.  2.3.  Decomposition  Site  on 
m-Fsce  of  Pure  AP  After  140  Min¬ 
utes  Decomposition  at  226°C  snd 
20  mm  Hg. 


The  decomposition  under  atmospheric  pressure  does  not  change  the 
general  appearance  of  the  decomposition  sites.  Figure  2.4  shows  part  of 
sn  m-fsce  sfter  90  minutes  decomposition  at  226°C  and  Fig.  2.5  psrt  of  a 
C-fsce  after  120  minutes  decomposition. 

However,  some  differences  from  vacuum  decomposition  csn  be  seen. 
First,  the  size  of  the  holes  is  smaller  although  the  crystals  decomposed 
about  twice  ss  long  as  under  vacuum.  This  is  probably  due  to  s  sup¬ 
pression  of  sublimstion  st  atmospheric  pressure.  Second,  there  is  a 
great  number  of  fine  crscks  not  visible  on  vscuun  decomposed  samples, 
snd  third,  there  sre  spots  of  elevated  surfsce.  These  spots  sre  concen¬ 
trated  at  the  rim  of  an  m-face  decompoaition  site  whereas  they  are  abun¬ 
dant  over  the  whole  streak  on  s  C-fsce. 

Aa  a  by-product  of  the  decomposition  at  both  20  mm  Hg  and  st  atmos¬ 
pheric  pressure  a  condensste  was  obtsined  on  a  cool  (about  50°C)  quartz 
plate.  Figure  2.6  shows  a  SEM  photograph  of  the  condensste  obtained 
during  the  decomposition  of  an  AP  crystal  at  one  atmosphere  of  air.  The 
condensste  wss  in  the  esrly  stsge  of  development.  An  x-ray  powder  pat¬ 
tern  showed  the  condensate  to  be  pure  AP.  This  indicates  that  there  waa 
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FIG.  2.4.  Decomposition  Site  FIG.  2.S.  Decomposition  Site 

on  m-Face  of  Pure  AP  After  90  on  C-Face  of  Pure  AP  After  120 

Mlnutea  Decomposition  st  226°C  Minutes  Decomposition  at  226°C 

and  Atmospheric  Pressure.  and  Atmospheric  Pressure. 


FIG.  2.6.  Condensate  of 
Pure  AP  on  Quartz  Disk. 
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an  appreciable  sublimation  of  AP  even  at  atmospheric  pressure.  No  con¬ 
densate  is  obtained  if  the  whole  reaction  vessel  (volume)  is  kept  at  de¬ 
composition  temperature.  In  that  case,  the  part  of  AP  which  sublimes 
subsequently  decomposes  in  the  gas  phase  or,  heterogeneously,  on  the  hot 
surface  (walls)  of  the  oven.  In  view  of  continuing  discussions  of  con¬ 
densed  and/or  gas  phase  reactions  during  the  decomposition  of  AP,  and  in 
light  of  the  above  observations,  an  effort  should  be  made  to  obtain  data 
simultaneously  for  reactions  occurring  in  the  condensed  phase  and  away 
from  the  condensed  phase  of  decomposing  AP. 

2.2.2.  Interfacial  Decomposition  of  Doped  AP  Crystals 

Work  in  this  area  as  reported  in  Ref.  1  has  been  expanded  to  in¬ 
clude  two  additional  concentrations  of  MnO^  in  the  doped  AP  single  crys¬ 
tals.  The  procedure  is  described  and  results  for  pure  AP  crystals  were 
reported  in  Ref.  1. 

Single  AP  crystals  containing  0.009  weight  percent,  0.04  weight 
percent,  and  0.11  weight  percent  MnO^"  were  decomposed  isothermally  under 
atmospheric  pressure.  The  advance  of  the  decomposition  interface  as  a 
function  of  time  was  obtained  by  cleaving  the  crystals  and  measuring  the 
thickness  of  the  residue  layers  at  different  stages  of  decomposition. 

The  total  weight  loss  after  a  complete  partial  decomposition  was  between 
35  and  40%  of  the  original  weight  as  compared  with  35%  for  pure  AP  and 
about  30%  measured  by  Bircumshaw  and  Newman  (Ref.  5). 

As  In  the  case  of  pure  AP,  it  was  found  that  the  decomposition 
interface  advanced  linearly  with  time  into  the  crystal.  Figures  2.7 
through  2.12  show  data  of  m-  and  c-face  decomposition.  The  decomposition 
parallel  to  C-faces  is  about  10-20%  slower  than  the  decomposition  paral¬ 
lel  to  w-faces.  With  the  rate  values  evaluated  from  Fig.  2.7  through 
2.12  one  can  calculate  the  fractional  decomposition  at  a  given  time  for 
crystals  which  were  cleaved  to  the  shape  of  a  parallelepiped.  Of  more 
interest,  however,  are  the  kinetic  parameters  of  the  isothermal  decom¬ 
position  which  can  be  obtained  from  Arrhenius  plots  of  rates  versus  1/T. 
Figure  2.13  shows  the  diagrams  for  m-face  decomposition  and  Fig.  2.14 
for  c-face  decomposition.  The  kinetic  parameters  calculated  from  these 
figures  are  listed  in  Table  2.1.  Values  reported  for  pure  AP  single 
crystals  (Ref.  1)  were  6.5  x  10^  cm/sec  for  the  preexponential  factor 
and  21.9  kcal/mole  for  the  activation  energy  (average  values  for  m-  and 
C-faces).  All  of  the  doped  AP  had  preexponencial  factors  and  activation 
energies  smaller  than  the  values  for  pure  AP.  Accordingly,  at  higher 
temperatures  and  under  the  assumption  that  the  decomposition  mechanism 
stays  the  same,  one  would  expect  smaller  burning  rates  for  doped  AP 
because  the  preexponential  term  becomes  more  and  more  dominant  and  the 
Arrhenius  plots  intersect  at  relatively  low  temperatures.  This  was  in 
fact  observed  (Section  2.3.2),  but  it  does  not  infer  that  the  combustion 
of  AP  obeys  the  same  or  similar  mechanisms  as  the  low  temperature  de¬ 
composition.  This  would  be  a  gross  simplification  of  a  very  complex 
process. 
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FIG.  2.7.  0.009  Wt.Z  Mn04. 

m-Face  Decompoaition. 


FIG.  2.9.  0.04  Wt.Z  MnO“. 

m-Face  Decomposition. 


FIG.  2.11.  0.11  Wt.Z  MnO^. 

m-Face  Decompoaition. 


FIG.  2.8.  0.009  Wt.Z  Mn04. 

C-Face  Decompoaition. 


FIG.  2.10.  0.04  Wt.Z  Mn04. 

C-Fsce  Decomposition. 


FIG.  2.12.  0.11  Wt.Z  MnO^. 

C-Face  Decomposition 


FIG.  2.7-2.12.  Penetration  of  the  Decomposition  Interface  Into  NH4Mn04 
Doped  Ammonium  Perchlorate  Single  Crystals,  d  ■  distance  from  crystal 
surface  to  Interface. 
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FIG.  2.13.  Arrhenius  Plot  for 
NH^MnO^ -Doped  AP.  m-Face  De¬ 
composition. 


FIG.  2.14.  Arrhenius  Plot  for 
NH^MnO^-Doped  AP.  C-Face  De¬ 
composition. 


TABLE  2.1. 


KINETIC  PARAMETERS  FOR  THE  DECOMPOSITION  OF  NH4Mn04  OOPED 
AMMONIUM  PERCHLORATE  SINGLE  CRYSTALS 


SAMPLE 

PREEXP.  FACTOR 
{cm/wcl 

m-fac*  c-fact 

ACTIVATION  ENERGY 

kcal/mok 
m-fici  c-facc 

AP 

+  0.009  WT.  % 
MnOj 

0.9 

62.0 

12.9 

17.3 

AP 

+  0.04  WT.  % 
MnO^ 

47.3 

39.0 

16.8 

16.7 

AP 

♦  0.11  WT.  % 
MnOj 

41800 

2.5 

21.4 

14.2 

A  theoretical  interpretation  of  the  results  has  not  been  tried  yet 
in  view  of  the  uncertainty  concerning  the  nucleation  rate  and  rate  of 
nuclei  growth.  In  preliminary  experiments  using  a  hot  stage  optical 
microscope  it  was  found  that  the  nucleation  rate  of  doped  AP  is  much 
higher  than  that  of  pure  AP.  In  addition,  the  number  of  nucleation 
sites  ia  orders  of  magnitude  greater  for  doped  than  for  undoped  AP  if 
both  are  compared  after  the  same  decomposition  time.  This  may  affect 
the  preexponentiai  factors  as  weil  as  the  activation  energies  and  a 
careful  atudy  of  nucleation  and  nuclei  growth,  both  for  pure  and  doped 
AP,  will  be  necessary  to  give  an  interpretation  of  the  variety  of  kine¬ 
tic  parameters  available. 
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2.3.  SINGLE  CRYSTAL  DEFLAGRATION  STUDIES 
2.3.1.  CombuBtlon  In  Reactive  Atmospheres 

In  earlier  work  (Ref.  1,  6  and  7)  evidence  was  reported  which  Indi¬ 
cated  that  a  thin  layer  of  melt  was  preaent  on  the  aurface  of  pure  AP 
cryatals  undergoing  aelf-auatalned  deflagration.  The  atudlea  have  been 
extended  to  determine  the  extent  of  occurrence  of  the  melt  and  lta  Im¬ 
portance  to  the  deflagration  procesa.  The  uae  of  reactive  atmoapherea 
auch  as  methane  and  ammonia  have  allowed  atudlea  to  be  made  at  reduced 
preBaurea  correapondlng  to  regreaalon  ratea  Intermediate  between  those 
of  self-deflagration  and  low  temperature  decompoaltlon. 

Single  crystals  of  pure  AP  and  crystals  doped  with  various  concen¬ 
trations  of  NH^MnO^  and  KCIO4  were  grown  in  the  apparatus  described  in 
Ref.  1.  Figure  2.15  shows  crystals  typical  of  those  used  in  the  present 
work.  Samples  were  burned  In  the  flushed  window  bomb  described  In 
Section  3.2.  High  speed  motion  pictures  were  taken  to  study  the  quali¬ 
tative  aspects  of  combustion  and  to  determine  the  burning  rates.  For 
photographic  purpoaea  pure  AP  and  K+-doped  cryatals  were  Illuminated  by 
a  2,500-watt  xenon  lamp.  The  purple  MnO^-doped  AP  crystals  were  photo¬ 
graphed  with  transmitted,  low  Intensity  polarized  light  to  prevent  heat¬ 
ing  the  crystals  by  optical  radiation  absorption. 


FIG.  2.15.  Single  Crystals 
Typical  of  Those  Used  In  the 
Present  Work.  At  left  Is  a 
MnO^-doped  crystal  which  has 
deep  purple  color  and  slightly 
modified  growth  habit.  At  right 
e  pure  AP  cryetal  snd  above 
amall,  cleaved  segments  of  Mn04~ 
doped  AP  In  various  concentra¬ 
tion  . 


High  speed  motion  pictures  of  the  deflagration  In  methane  at  pres¬ 
sures  between  etmospheric  and  100  psla  were  of  perticular  Interest. 

Under  these  conditions  a  variety  of  surface  heating  conditions  end, 
hence,  regression  retes  was  possible.  Figure  2.16  shows  s  pure  AP 
crystal  burning  st  15  psle  In  methane.  At  this  pressure  the  surface 
regressed  with  a  pronounced  coning  due  to  flame  spread  down  the  sldee 
of  the  crystal;  the  upper  portion  of  the  coned  eurface  had  a  regression 
rate  of  approximately  0.03  ln/sec.  The  surface  had  e  froety,  translu¬ 
cent  eppearsnce.  Quenched  samples  st  these  lower  pressures  were  obtelned 
by  rspidly  flooding  the  combustion  bomb  with  nitrogen.  Figure  2.17  Is  a 
scanning  electron  microscope  (SEM)  photograph  of  the  surface  of  a 
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FIG.  2.16.  Single  Frame  of 
Motion  Picture  of  Pure  AP 
Burning  In  Methane  Atmo¬ 
sphere  at  15  psla.  Crystal 
approximately  0.5  cm  wide. 


FIG.  2.17.  Scanning  Elec¬ 
tron  Microscope  Photomicro¬ 
graph  of  the  Burned  Surface 
of  Pure  AP  Quenched  from 
Combustion  In  15  psla 
methane. 


quenched  sample  that  was  burned  in  this  intermediate  regression  rate 
domain.  The  surface  structure  that  results  from  this  mode  of  combustion 
gives  the  appearance  of  a  crystal  structure  but  differs  from  the  struc¬ 
ture  of  an  "undisturbed"  AP  crystal.  Figure  2.18  shows  a  pure  AP  crys¬ 
tal  burning  at  65  psla.  At  this  pressure  the  surface  was  clear  and 
transparent  and  resembled  melting  ice.  Occasional  boiling  activity 
could  be  clearly  seen  on  the  surface.  At  100  psia  in  methane  movies 
taken  at  a  magnification  of  AX  of  AP  crystals  approximately  250  micron 
thick  showed  a  definite  "flowing"  motion  associated  with  portions  of 
the  burning  surface  as  well  as  patches  of  rapidly  boiling  "froth".  It 
was  not  possible  to  quench  samples  under  these  conditions  of  higher 
pressure. 


11 


FIG.  2.18.  Single  Frame  of 
Motion  Picture  of  Pure  AP  Burn¬ 
ing  In  Methane  Atmosphere  st 
65  psla.  Crystal  approximately 
0.5  cm  wide. 


The  K^-doped  AP,  as  observed  in  the  movies,  behaved  eaaentlally  like 
pure  AP  when  burning  in  methane  except  that  the  tendency  for  surface 
boiling  waa  more  pronounced  and  eaally  observed. 


In  low  pressure  (15  to  50  paia)  methaoe  the  KnO^-doped  crystals  were 
observed  to  burn  with  a  white  porous  material  covering  the  regressing 
surface.  Figure  2.19  is  a  SEM  photograph  of  the  aurface  of  a  MnO^-doped 
cryatal  quenched  at  low  preaaure.  This  structure  is  very  similar  to  the 
water  vapor  "rejuvinated"  low  temperature  decomposition  residue  shown  In 
Fig.  2.20.  At  higher  pressure  (200  psla)  the  poreua  material  appeared 
to  undergo  a  cyclic  accumulation-rapid  consumption  process  which  gave  a 
strsnge  pulsating  character  to  the  deflagration  process. 


FIG.  2.19.  Scanning  Electron  Mi¬ 
croscope  Photomicrograph  of  the 
Burned  Surface  of  MnO^-Doped  AP 
Quenched  from  Combustion  In  15 
psla  Methane. 


FIG.  2.20.  Scanning  Electron  Mi¬ 
croscope  Photomicrograph  of  the 
Residue  of  Partially  Decomposed 
Pure  AP  (st  220°C)  Which  has  Been 
Rejuvinated  by  Exposure  to  Water 
Vapor. 
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Low  speed  motion  pictures  (24  to  100  fps)  were  taken  utilizing  the 
self-illumination  of  the  AP  decomposition  flame.  At  atmospheric  pres¬ 
sure  in  methane  a  pink-orange  flame  approximately  175  microns  thick  was 
observed  (by  visual  observation  the  flame  color  was  the  same  as  was  ob¬ 
served  with  single  crystals  burning  in  nitrogen  at  elevated  pressure). 
The  decomposition  flame  appeared  to  commence  right  at  the  condensed 
phase  surface. 

In  addition  to  the  samples  burned  in  methane,  several  pure  AP 
crystals  were  burned  in  an  atmosphere  of  ammonia  at  atmospheric  pres¬ 
sure.  The  combustion  behavior  was  essentially  the  same  as  in  methane 
atmospheres.  It  was  not  possible  to  burn  samples  at  elevated  pressure 
in  ammonia  since  the  crystals  rapidly  dissolved  in  the  ammonia. 

2.3.2.  Deflagration  at  High  Pressure 

Additional  measurements  were  made  of  the  deflagration  rates  for 
single  crystals.  The  rates  were  meaaured  from  high  speed  motion  pic¬ 
tures  of  the  combustion  and  are  shown  in  Fig.  2.21.  The  previously 
reported  AP  burning  rate  data  (Ref.  6)  were  repeated  and  verified  with 
larger  cross-section  (0.5  cm  on  a  side)  crystals. 


FIG.  2.21.  Deflagration  Rate  of 
Pure  and  Doped  AP  Single  Crystals. 
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The  burning  rates  for  the  single  AP  crystsls  doped  with  K+  and 
three  concentrations  of  MnO^  exhibited  one  consistent  feature  even 
though  there  was  Inconsistency  in  some  of  the  rate-concentration  dats; 
that  feature  was  that  the  burning  rates  for  the  doped  crystals  were  al¬ 
ways  below  the  pure  AP  rate  over  the  pressure  range  tested.  In  addition, 
the  K+-doped  AP  had  a  pressure  exponent  of  approximately  1.0  as  compared 
to  approximately  0.77  for  pure  AP.  Motion  pictures  of  this  material 
undergoing  deflagration  showed  definite  melting  snd  boiling  on  the  sur¬ 
face  with  occasional  ejection  of  molten  material  (presumed  KC1)  from  the 
surface. 

In  the  case  of  the  MnOJ'-doped  crystals,  the  burning  rate  data  at  the 
lowest  concentrations  was  too  scattered  to  warrant  drawing  conclusions 
concerning  rate-concentration  trends.  In  addition,  it  was  difficult  to 
precisely  measure  the  MnO^  concentration.  The  highest  (0.04X)  concentra¬ 
tion  Mn04-doped  AP  which  would  undergo  self-deflsgration  displayed  a  much 
lower  burning  rate  so  it  is  presumed  that  increasing  MnO^  concentration 
results  in  lowering  the  single  crystal  deflagration  rate.  Light  trans¬ 
mission  through  the  gaseous  combustion  products  was  observed  to  be  much 
higher  for  the  MnOJ-doped  AP  than  for  the  pure  AP.  Figure  2.22  shows  a 
comparison  of  doped  versus  pure  AP  photographed  under  the  same  test  con¬ 
ditions  with  transmitted  Illumination  by  polarized  light. 


(a)  (b) 

FIG.  2.22.  Single  Frame  of  Motion  Picture  of  (s)  Pure  AP  and 
(b)  AP  Doped  With  0.006X  MnCJ".  Both  samples  deflagrating  st 
800  psia,  photographed  with  transmitted,  polarized  illumination. 
Note  difference  In  optical  density  of  combustion  products. 
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The  low  pressure  deflagration  limit  for  all  of  the  doped  crystals 
was  much  higher  than  that  of  the  pure  material  (275  psia.  Ref.  6)  with 
MnO^T  being  apparently  more  effective  in  raising  the  limit  than  K+. 
Friedman,  et  al  (Ref.  9),  also  observed  that  NaMn04  .  3H2O)  and  Mnt^ 
raised  the  deflagration  limit  of  pressed  AP  pellets.  The  low  pressure 
deflagration  limits  in  the  present  work  were  found  to  be  approximately 
600  psia  for  the  0.5%  KClO^-doped  AP  and  800  psia  for  the  0.003  to 
0.006%  MnO^-doped  AP.  The  AP  doped  with  0.04%  MnO^  could  not  be  made  to 
deflagrate  below  about  2000  psia. 

2.3.3.  Interpretation  and  Discussion 

The  work  here  and  in  Ref.  1  has  revealed  several  qualitative  as¬ 
pects  of  combustion  that  were  previously  unknown,  or  were  matters  of 
speculation.  The  most  important  aspect  is  the  decisive  evidence  for  a 
frothy  liquid  on  the  surface  of  the  burning  AP.  The  exact  chemical 
nature  of  the  liquid  is  not  known,  although  x-ray  diffraction  tests  on 
residues  indicate  only  AP.  It  is  not  surprising  that  a  liquid  layer 
shauld  be  present  at  the  high  surface  heating  rates  associated  with  self¬ 
deflagration  of  crystalline  AP.  The  temperature  gradient  resulting  from 
this  heat  flux  is  of  such  magnitude  that  it.  could  raise  the  temperature 
of  the  subsurface  AP  to  its  melting  point  in  less  time  than  it  would 
take  to  gasify  the  material.  Thus,  a  melting  point  would  be  encountered 
which  would  not  be  observed  at  lower  heating  rates,  i.e.,  such  as 
associated  with  differential  thermal  analysis,  where  the  sample  would 
gasify  before  the  melting  point  was  reached. 

The  results  suggest  that  the  melt  plays  an  important  role  in  self¬ 
deflagration  of  the  AP.  There  has  been  continuing  difficulty  in  recon¬ 
ciling  the  amount  of  heat  required  for  the  initial  decomposition  of  AP 
with  the  amount  of  heat  that  could  be  supplied  by  the  AP  flame,  because 
the  heat  balance  dictates  a  flame  standing  only  a  few  mean  free  paths 
from  the  surface  (assuming  a  one-dimensional  regressing  surface).  On 
the  other  hand,  a  frothy  melt  on  the  surface  provides  opportunity  for 
both  heterogeneous  and  gas  phase  reactions  in  the  froth,  with  a  very 
large  area  for  heat  transfer.  Thus,  it  seems  plausible  to  assume  that 
the  liquid  froth  is  an  essential  feature  of  the  self-deflagration  of  AP, 
and  that  the  low  pressure  deflagration  limit  corresponds  to  attainment 
of  heating  rates  so  low  that  surface  material  gasifies  without  reaching 
a  melting  temperature.  Under  these  conditions,  a  major  portion  of  the 
exothermic  reactions  would  occur  only  in  the  gas  phase  flame,  which 
apparently  does  not  have  any  equilibrium  position  satisfying  heat  balance 
requirements . 

Microscopic  examination  of  quenched  samples  indicate  no  optically 
detectable  changes  in  the  AP  below  the  melt  region  other  than  the  phase 
transition  region.  This  implies  that  there  is  probably  no  significant 
contribution  to  the  heat  release  during  deflagration  by  solid  phase  re¬ 
actions.  This  is  not  really  surprising  because  the  decomposition  re¬ 
actions  of  solids  include  an  induction  period  followed  by  a  growth  and 


15 


NWC  TP  4478 


decay  period.  The  induction  period  may  be  long  or  at  least  comparable 
to  the  time  it  takes  for  the  very  thin  thermal  wave  associated  with  com¬ 
bustion  to  traverse  a  given  layer  of  material.  Ammonium  perchlorate 
does  decompose  in  a  combustion  wave,  and  on  a  molecular  scale  the  actual 
decomposition  mechanism  may  be  the  same  as  for  the  "low  temperature" 
thermal  decomposition,  e.g.,  a  proton  transfer.  The  kinetics  of  the 
overall  process  are  expected  to  be  different  if  they  occur  in  s  melt  as 
opposed  to  a  crystalline  solid. 

It  aeems  reasonable  to  asaume  then,  that  slteration  of  the  thermal 
stability  of  the  surface  melt  or  auppresaing  its  formation  would  produce 
a  corresponding  change  in  the  deflagration  behavior  of  AP.  Earlier  work 
(Ref.  9)  showed  that  a  rapidly  heated  mixture  of  granular  AP  and  lithium 
fluoride  produced  a  residue  mixture  of  AP-lithium  perchlorate.  Markowitz 
and  Harris  (Ref.  10)  have  shown  that  the  decomposition  rste  of  AP  in 
solution  with  lithium  perchlorate  is  lower  than  that  of  pure  AP  at  corre¬ 
sponding  temperatures.  It  is  possible  that  a  similar  process  may  occur 
on  the  surface  of  deflagrating  crystals  doped  with  the  K+  ion  resulting 
in  the  formation  of  an  AP-potaasium  perchlorate  solution.  There  is  also 
an  accumulation  of  liquid  balls  (probably  KC1,  a  stable  product  of  KCI4 
decomposition)  on  the  burning  surface  of  K+  doped  crystals  which  may 
affect  the  deflagration  rate.  In  the  case  of  MnO^-doped  crystals  both 
the  deflagration  limit  and  the  burning  rate  were  Influenced  tremendously 
by  small  concentrations  of  this  ion.  The  preasure  deflagration  limit 
is  Increased  accordingly  along  with  a  reduced  burning  rate. 


2.4.  SUMMARY 

The  present  work  (including  Ref.  1)  has  revealed  rather  decisively 
that  the  physical  character  of  the  surface  residue  of  AP  is  significantly 
dependent  on  the  imposed  heating  rate.  Under  conditions  of  isothermal 
decomposition  the  AP  apparently  sublimes  and  there  is  no  evidence  of  a 
surface  melt.  Under  conditions  of  rapid  heating,  a  surface  melt  is  ob¬ 
served  and  in  the  case  of  self-deflagration  the  surface  melt  is  frothy 
and  may  contribute  significantly  to  the  heat  release  necessary  to  sustain 
deflagration.  Also,  the  substitution  of  foreign  iona  in  the  AP  crystal 
structure  in  amall  quantities  reduced  the  activation  energy  of  the  iso¬ 
thermal  decomposition  snd  also  altered  the  combustion  behavior  of  the 
crystal  either  by  altering  the  thermal  stability  of  the  aurface  melt 
and/or  the  gas  phase  kinetics. 
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3.  BEHAVIOR  OF  ALUMINUM  IN  COMPOSITE  PROPELLANTS 


Certain  aspects  of  aluminum  agglomeration  in  the  combustion  zone  of 
composite  propellants  were  described  in  detail  in  the  previous  progress 
report  (Ref.  1).  During  this  reporting  period,  emphasis  has  been  placed 
on  investigating  the  mechanism  of  agglomeration.  In  particular,  the 
photographic  technique  has  been  improved  to  increase  the  resolution  of 
the  high  speed  motion  pictures,  and  quenched  propellant  samples  have  been 
examined  in  detail  with  a  scanning  electron  microscope  (SEM)  to  investi¬ 
gate  the  early  (preignition)  phases  of  aluminum  agglomeration. 


3.1.  EXPERIMENTAL  TECHNIQUE  DEVELOPMENTS 

3.1.1.  Clnephotomlcrographv  Improvements 

The  basic  arrangement  for  high  speed  motion  picture  photography  with 
magnification  has  been  described  in  Ref.  1.  However,  improvements  have 
been  made  in  each  of  the  three  major  components  of  the  experimental  setup. 

The  improved  arrangement  is  shown  in  Fig.  3.1.  The  original  Fastax 
camera  was  replaced  by  a  Hycam  16  mm  high  speed  motion  picture  camera, 
Model  K3004E.  This  camera  provides  more  uniform  exposure  over  the  entire 
film  frame,  more  versatility  in  speed  selection  and,  due  to  its  Internal 
optics,  provides  a  somewhat  larger  magnification  than  the  Fastax.  The 
original  xenon  arc  external  illumination  arrangement  was  replaced  with  a 
Strong  X-16  xenon  arc  assembly  which  includes  a  2,500  watt  xenon  arc 
lamp,  dichroic  main  mirror  and  external  focusing  controls.  This  assembly 
permits  better  control  of  external  illumination  and  can  supply  consider¬ 
ably  more  intensity,  if  desired,  than  the  previous  arrangement.  The 
greatest  improvement  has  been  made  in  the  window  bomb  itself.  The  new 
bomb  provides  complete  flushing  of  the  interior  without  recirculation, 
permits  more  versatility  in  the  mounting  of  propellant  samples  (Fig.  3.2), 
and  provides  better  mechanical  support  for  the  fused  quartz  windows. 

Figure  3.3  illustrates  the  details  of  combustion  that  can  now  he 
resolved  with  the  improvements  in  the  experimental  arrangement.  In  the 
first  picture,  an  accumulation  of  several  hundred  15y  aluminum  particles 
have  melted  together  to  form  a  sphere  of  molten  aluminum,  as  indicated 
by  the  highlight  on  die  sphere.  Some  unmelted  material,  including 
several  original  aluminum  particles,  is  visible  on  the  left  side  of  the 
agglomerate.  In  the  second  picture  (one  millisecond  later)  most  of  the 
remaining  material  on  the  agglomerate  has  melted  except  for  the  crusty. 
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FIG.  3.3  Selected  liamru  from  hifth  a|)PfJ  motion  picturr  of  com- 
ponitr  propc'IUnl  containing  18%  aluminum  burning  at  100  pni. 
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irregular  appearing  aluminum  oxide  surface.  The  aluminum  oxide  "smoke" 
indicates  the  developing  reaction  region  of  the  aluminum  combustion.  At 
five  milliseconds  the  droplet  shows  a  well-developed  combustion  envelope 
which  flows  outward  with  the  propellant  combustion  products.  In  the 
final  picture  the  spherical  droplet  appears  to  be  burning  vigorously  and 
highlights  from  the  liquid  surface  are  visible  through  the  combustion 
envelope.  The  disc  centered  on  the  sphere  in  this  picture  is  a  luminous 
cap  which  probably  contains  most  of  the  preignition  oxide  of  the  original 
aluminum  particles.  Such  oxide  caps  are  readily  discernible  in  the 
photographs  in  Ref.  1  (Fig.  3.3).  The  final  picture  also  shows  the  het¬ 
erogeneity  of  the  propellant  surface  underneath  the  aluminum  agglomerate. 
The  motion  pictures  show  this  area  to  be  almost  entirely  "uncluttered" 
with  pyrolysis  residue.  Assuming  that  the  residue  is  a  thermal  insulator, 
energy  from  the  gas  phase  could  be  more  easily  transferred  to  this  '"clean" 
spot  and  could  result  in  a  momentary  increase  in  energy  release  in  that 
localized  area.  This  possibility  along  with  the  energy  release  following 
ignition  of  the  aluminum  agglomerate  quite  conceivably  could  play  an 
important  role  in  driving  gas  oscillations. 

3.1.2.  Propellant  Quench  Technique 

Quenching  propellants  by  rapid  depressurization  of  the  window  bomb 
was  performed  earlier  in  this  program  (Ref.  11).  High  speed  motion 
pictures  taken  during  the  quench  indicated  that  the  majority  of  aluminum 
in  the  later  stage  of  agglomeration  was  removed  from  the  surface  during 
the  rapid  depressurization.  In  order  to  preserve  the  large  accumulates 
during  quench  for  later  microscopic  study,  a  new  quench  technique  was 
employed.  A  freshly  cut  piece  of  propellant  (approximately  a  1/4  in.  cube 
but  with  the  top  surface  cut  at  a  30°  angle  to  the  horizontal)  was  placed 
on  a  clean  metal  plate  and  ignited  by  a  hot  Michrome  wire  in  the  pressur¬ 
ized  bomb.  As  the  burning  surface  approached  the  metal  plate,  the  com¬ 
bustion  was  thermally  quenched.  A  small  portion  of  unburned  propellant 
as  well  as  aluminum  particles  in  various  stages  of  agglomeration  remained 
on  the  plate.  Plates  of  copper,  aluminum,  and  stainless  steel  were  eval¬ 
uated  as  heat  sinks;  the  stainless  steel  plates  were  most  satisfactory 
in  that  a  larger  portion  of  unburned  propellant  remained  after  quench. 

Some  residuums  have  been  examined  with  a  scanning  electron  microscope 
(SEM)  and  the  results  are  discussed  in  later  sections. 


3.2.  DETAILED  ASPECTS  OF  ALUMINUM  AGGLOMERATION 

The  preignition  behavior  of  the  aluminum  in  composite  propellants 
has  received  increased  attention  during  this  reporting  period  and  the 
more  important  processes  are  reasonably  well  understood.  An  exploratory 
experiment  which  yielded  the  greatest  insight  into  the  agglomeration 
mechanism  was  performed  by  dropping  lOOp  spherical  aluminum  particles 
(in  air)  onto  a  sapphire  disk  heated  by  a  torch  to  approximately  1,400°C 
(Ref.  12)  and  subsequently  examining  the  particles  using  a  scanning 
electron  microscope  (SEM).  A  series  of  resulting  SEM  photographs  with 
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increasing  magnification  are  shown  in  Fig.  3.4.  The  important  feature 
shown  in  Fig.  3.4a  is  that  two  aluminum  particles  are  rigidly  connected 
by  a  "bridge"  of  material.  This  interaction  of  the  particles  occurred 
at  an  environmental  temperature  (<  1,400  C)  well  below  the  melting  point 
of  aluminum  oxide  (2,050°C),  even  though  the  original  aluminum  particles 
had  an  oxide  coating.  Examination  at  higher  magnifications  shows  numer¬ 
ous  ridges  approximately  l/2p  in  width  which  resemble  welding  joints. 

These  observations  on  the  microscopic  scale  probably  explain  the  details 
of  aluminum  agglomeration  in  propellants  and  are  interpreted  in  the 
following  manner.  As  the  aluminum  particles  are  heated  (in  this  experi¬ 
ment  by  conduction  from  the  1,400°C  sapphire  disk),  the  stresses  produced 
by  the  differential  thermal  expansion  of  the  aluminum  and  aluminum  oxide 
caused  the  aluminum  oxide  covering  to  fracture.  The  particles  in  Fig.  3.4 
were  evidently  heated  to  above  the  melting  point  of  aluminum  and  hence, 
liquid  aluminum  started  to  flow  through  the  cracks.  The  aluminum  was 
rapidly  oxidized,  thus  producing  the  "weld  joint"  appearance.  The  heat 
liberated  by  the  aluminum  oxidation  conceivably  could  produce  additional 
stresses  and  cracking  if  radiative  heat  losses,  for  instance,  were  not 
too  great.  As  in  Fig.  3.4,  where  the  particles  are  in  contact,  a  "bridge" 
can  be  formed  by  the  cracking  of  oxide  and  flow  of  molten  aluminum.  In 
Fig.  3.4d  it  appears  that  additional  "thermal  cracking"  occurred  after 
the  bridge  was  formed.  This  is  indicated  by  the  two  ridges  that  run 
across  the  bridge. 

It  is  felt  that  this  phenomena  is  the  first  step  in  the  agglomer¬ 
ation  (or  sintering  together)  process.  Combustion  residues  obtained 
from  thermal  quenching  (Section  3.1.2)  were  examined  with  the  SEM  for 
evidence  of  the  "bridging”  phenomena.  Figure  3.5a  is  a  series  of  SEM 
photographs  of  an  agglomerate  from  a  composite  propellant  containing  16% 
of  15u  diameter  aluminum.  The  successively  larger  magnifications  show 
the  pattern  of  the  aluminum  flowing  to  form  a  large  agglomerate.  For 
comparison,  Fig.  3.5b  shows  an  analogous  situation  of  agglomerate  for¬ 
mation  with  1 OOp  aluminum  particles  rapidly  heated  to  1,400  C  on  a 
sapphire  disk  in  atmospheric  air.  The  details  of  oxide  cracking  and 
healing  are  not  at  all  evident  in  the  agglomerate  from  the  propellant. 

This  agglomerate  was  formed  in  a  more  complicated  environment  and  un¬ 
doubtedly  there  are  numerous  pyrolysis  and  oxidation  products  present 
which  could  greatly  modify  the  appearance  of  the  agglomerate.  However, 
efforts  are  continuing  to  conclusively  demonstrate  this  "bridging" 
phenomena  because  of  the  importance  of  relating  the  metal  particle  be¬ 
havior  with  behavior  in  propellants. 


3.3.  INTERPRETATION  AND  SUMMARY 

A  detailed,  qualitative  description  of  the  behavior  of  aluminum  in 
composite  propellants  can  now  be  made  based  on  (1)  high  speed  motion 
pictures  of  aluminized  propellants,  especially  the  observations  of  a 
liquid  phase  on  the  propellant  surface  for  the  early  stages  of  accumu¬ 
lation,  (2)  on  interpretation  of  and  deductions  from  the  observations  of 
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FIG.  3.4.  SEM  Photographs  of  "Bridgea”  Joining  lOOp 
Aluminum  Particles.  Aluminum  particles  were  dropped 
onto  sapphire  plate  heated  to  1I400°C  In  air. 

the  behavior  of  rapidly  heated  single  particles  in  Section  3.2  and 
similar,  though  less  conclusive,  studies  of  residual  accumulates  from 
propellant  quench  tests,  and  (3)  on  high  speed  motion  pictures  (Section 
3.1.1  and  Ref.  1  and  13)  of  burning  propellants  as  exemplified  by  the 
four  frames  shown  In  Fig.  3.3. 

The  sequence  of  events  in  the  agglomeration  of  aluminum  is  believed 
to  be  as  follows.  The  original  aluminum  particles  with  thin  oxide  coat- 
lngs  are  distributed  In  the  binder  of  the  propellant.  As  the  combustion 
front  movea  through  the  propellant,  the  temperature  of  the  aluminum  par¬ 
ticles  begins  to  rlae  as  the  thermal  wave  proceeds.  As  the  binder  and 
the  ammonium  perchlorate  pyrollze,  the  aluminum  recedes  with  the  surface 
until  it  comes  into  contact  with  other  aluminum  particles.  The  decom¬ 
posing  surface  of  the  propellant  is  liquid  and  thus  the  Individual 
aluminum  particles  are  readily  retained  on  the  surface.  As  the  aluminum 
particles  accumulate  in  the  melt,  they  are  heated  by  conduction  to  at 
least  the  local  surface  temperature  (typically  on  the  order  of  400°C- 
700°C)  and  probably  reach  a  somewhat  higher  temperature  due  to  oxidative 
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self-heating.  In  this  situation  with  the  aluminum  particles  in  contact 
with  one  another  and  with  the  aluminum  at  a  temperature  above  its  melting 
point,  the  adjacent  particles  interact  as  described  in  Section  3.2.  The 
molten  aluminum  flows  through  the  thermal  expansion  cracks  in  the  aluminum 
oxide  ahell  and  is  probably  quickly  oxidized.  When  cracks  roughly  coin¬ 
cide  with  the  contact  points  with  adjacent  particles,  the  flowing  aluminum 
can  "wet"  the  adjacent  particle.  Thua,  the  aluminum  particles  are  joined 
together  at  their  contact  points  by  aluminum  and  aluminum  oxide  "bridges". 
The  agglomerate  as  s  whole  can  be  described  as  numerous  drops  of  liquid 
aluminum  encapsulated  in  aluminum  oxide  shells.  These  drops  and  ahella 
are  linked  together  by  oxide  "bridges".  Figure  3.4  representa  this  type 
of  structure.  As  the  binder  and  AP  recede  beyond  this  region  of  the 
aluminum  agglomerate,  the  aluminum-aluminum  oxide  structure  haa  enough 
strength  to  physically  withstand  the  flow  of  binder  and  oxidizer  gases 
away  from  the  aurface  and,  nence,  remains  on  the  surface  until  it  ignites 
or  is  completely  undermined  and  drifts  away  with  the  flow  of  gsa  from  the 
aurface. 

In  the  case  where  the  agglomerate  is  not  undermined  and  removed, 
oxidative  self-heating,  oxide  cracking  and  aluminum  flow  continues  until 
the  molten  aluminum  from  one  ahell  joins  with  aluminum  from  another  and 
a  large  molten  ball  of  aluminum  has  formed.  The  latter  stage  of  the 
formation  of  a  large  molten  ball  of  aluminum  is  shown  in  the  "t"  portion 
of  Fig.  3.3.  The  next  picture  (t  +  1  msec)  ahows  the  ball  increasing  in 
aize  to  engulf  the  entire  agglomerate  contained  in  that  particular  binder 
"pocket".  These  two  pictures  also  show  that  the  aluminum  oxide  on  the 
large  molten  ball  ia  irregular  in  ahape  indicating  that  the  oxide  has 
not  yet  melted.  In  the  t  +  5  msec  picture  the  oxide  has  melted  (no  aharp 
irregularities)  and  vigorous  combustion  of  the  aluminum  has  begun  as 
indicated  by  the  pronounced  aluminum  oxide  smoke  that  has  developed.  The 
molten  aluminum  oxide  cap  characteristic  of  burning  aluminum  droplets  ia 
distinguishable  in  the  t  +  8  msec  picture. 

Future  effort  in  this  area  will  be  directed  at  obtaining  even  more 
conclusive  evidence  on  the  aluminum  sintering  process  and  the  transition 
to  a  large  molten  ball  in  experiments  closely  simulating  the  propellant 
combustion  situation. 
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4.  NONACOUSTIC  INSTABILITY 


4.1.  INTRODUCTION 

A  general  review  of  the  statua  of  nonacoustic  combustion  instabil¬ 
ity  (NAI)  in  aolid  propellants  was  presented  in  the  previous  report 
(Ref.  1).  The  principal  effort  during  this  reporting  period  has  been 
(i)  to  compare  nonacoustic  experimental  results  with  one-dimensional 
theories  and  (2)  to  examine  the  effect  of  oxidizer  particle  size  on  NAI. 

In  the  past,  theoretical  models  describing  the  interaction  of  the 
pressure  perturbations  with  the  combustion  zone  have  mainly  been  applied 
to  acoustic  instability  (for  example.  Ref.  14-17).  However,  the  combus¬ 
tion  instability  models  that  describe  the  dynamics  of  the  combustion 
zone  and  are  applicable  in  the  low  frequency  regime  should  apply  to  non¬ 
acoustic  as  well  as  acoustic  instability.  Several  recent  papers  (Ref. 
18-20)  have  discussed  NAI  in  terms  of  the  response  function  normally 
associated  with  acoustic  instability  which  evolved  from  the  above  models. 
In  the  previous  report  (Ref.  1)  the  analysis  relating  to  the  combustor 
was  separated  from  the  analysis  of  the  combustion  so  that  theoretical 
combustion  modela  could  be  combined  with  the  combustor  analysis  quite 
readily.  A  method  was  also  presented  for  reducing  NAI  data  (obtained 
from  L*  burners)  in  terms  of  the  response  function  and  the  phase  of  the 
combustion.  In  the  present  report,  some  of  the  same  NAI  data  are  com¬ 
pared  with  response  function  curves  calculated  from  models  which  are 
directly  applicable  to  the  low  frequency  regime  but  which  have  normally 
been  applied  to  acoustic  instability.  For  this  purpose  the  model  of 
Denison  and  Baum  (Ref.  16)  and  &  recent  model  proposed  by  Hart,  Farrell 
and  Cantrell  (Ref.  17)  have  been  chosen  as  being  representative  and  have 
been  utilized  to  compare  the  relative  applicability  of  these  types  of 
models . 

Experimental  experience  plus  the  understanding  obtained  from  studies 
such  as  outlined  above  have  advanced  the  level  of  understanding  of  NAI  so 
that  predictions  can  be  made  as  to  the  general  behavior  of  a  given  pro¬ 
pellant.  However,  there  are  many  instances  where  the  behavior  has  been 
different  than  expected.  For  example,  the  one-dimensional  models  predict 
that  for  a  given  propellant  and  a  given  set  of  conditions,  pressure  os¬ 
cillations  should  occur  at  one  particular  frequency.  However,  two  pre¬ 
vious  instances  are  knovn  (Ref.  I  and  21)  where  propellants  have  oscil¬ 
lated  at  two  different  frequencies,  under  essentially  the  same  conditions, 
and  in  each  case  one  of  the  frequencies  was  much  different  than  expected. 

A  close  scrutiny  of  the  propellant  characteristics  and  the  conditions 


2b 


NWC  TP  4478 


under  which  the  anomalies  occurred  has  led  to  the  belief  that  the  oxidizer 
particle  size  was  a  major  factor  producing  the  unexpected  frequency.  This 
provided  the  motivation  which  led  to  the  atudy  which  la  outlined  In 
Section  4.5. 


4.2.  NONACOUSTIC  COMBUSTOR  ANALYSIS 

In  NAI  It  appears  that  the  combustion  couplea  with  the  flow  out  of 
the  rocket  chamber  in  order  to  Induce  and  sustain  instability.  Thua,  the 
maas  balance  for  a  rocket  motor  containing  burning  propellant  and  ex¬ 
hausting  through  a  aonlc  nozzle  can  be  written 


pr  Afc  - 
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(4.1) 


where  c*  la  the  characteristic  velocity  for  the  products  of  combustion. 
Assuming  Ideal  gaa  conditions,  a  thin  combustion  zone,  s  spacewlse  Iso¬ 
thermal  gas,  and  the  density  of  the  propellant  to  be  much  greater  than 
that  of  the  gas,  Eq.  4.1  can  be  written  (For  a  more  complete  derivation 
of  Eq.  4.1  and  4.2,  aee  Ref.  18  and  22) 


(4.2) 


Applying  standard  perturbstion  techniques  and  considering  an  exponen¬ 
tially  growing,  sinusoidal  perturbation  such  as 
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Equation  4.2  reduces  to 


“l+OT.+idlT,  (4.4) 
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where  a  is  a  growth  coefficient  and  tcj<  represents  the  residence  time  of 
the  chamber  and  is  directly  proportional  to  L*.  The  left  hand  side  of 
Eq.  4.4  Is  the  response  of  the  mass  flow  rate  In  the  combustor  to  the 
pressure  of  the  combustor  and  is  equal  both  In  magnitude  and  phase  to  a 
similarly  defined  response  of  the  combustion.  All  of  the  terms  on  the 
right  hand  aide  of  Eq.  4.4  are  quantities  that  can  be  determined  ex- 
periemntslly  In  an  L*-type  of  burner.  Thua,  the  nonacoustic  response 
of  the  combustor  can  be  determined  experimentally  snd  compared  with  the 
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response  of  the  combustion  as  determined  from  models  such  as  those  of 
Denison  and  Baum  and  Hart,  Farrell  and  Cantrell. 

Figure  4.1  is  a  schematic  of  an  L*  burner  (Ref.  18,  21,  and  22) 
from  which  nonacoustic  instability  data  can  be  obtained.  Data  have  been 
obtained  for  two  propellants  and  reduced  to  the  format  suggested  by 
Eq.  4.4  in  order  to  make  a  comparison  with  one-dimensional  models.  Both 
are  composite  propellants  but  vary  widely  in  compositions  and  properties. 
A-35  is  a  polyurethane-AP  propellant,  containing  752  AP,  while  Utah-TF 
is  a  PBAA-AP-A1  propellant  containing  5%  Al,  182  PBAA,  22  copper  chromite 
burning  rate  catalyst  and  the  remainder  AP  (see  Ref.  1  and  21). 


MOPEU.ANT 


4.3.  COMBUSTION  DYNAMICS  MODELS 

The  model  of  Denison  and  Baum  (denoted  DB)  was  chosen  to  provide  a 
basis  of  comparison  for  several  reasons.  First  of  all  it  applies  specifi¬ 
cally  to  the  low  frequency  range  that  is  of  interest  in  NAI.  Secondly, 
as  has  been  pointed  out  by  Culick  (see  Appendix  A  of  Ref.  1),  the  form 
of  the  derived  response  function  is  extremely  simple,  containing  oniy  two 
parameters  (A  and  Bqjj)  and  the  independent  variable  D,  Culick  has  also 
shown  that  virtually  all  of  the  one-dimensionai  models  reduce  to  the  same 
function  as  the  DB  model.  The  model  of  Hart,  Farrell,  and  Cantrell 
(denoted  1IFC)  was  also  chosen  because  it  too  is  directly  applicable  to 


27 


NWC  TP  4478 


the  low  frequencies  of  interest  in  the  present  study,  and  because  it 
contains  the  effects  of  a  pressure  dependent  surface  reaction. 

4.3.1.  Comparing  the  Models 

The  assumptions  involved  in  both  of  the  analyses  are  very  similar. 
Both  assume  a  homogeneous,  one-dimensional  solid,  quasi-steady  gas  phase, 
and  are  applicable  in  the  same  low  frequency  range.  The  principal  dif¬ 
ference  in  the  models  (other  thsn  the  actual  method  of  analysis)  is  that 
Denison  and  Baum  arrive  at  a  burning  rate  lsw  by  starting  with  the  con¬ 
servation  equations  and  reducing  the  complexity  of  these  equations  with 
several  assumptions.  In  so  doing,  the  effects  of  diffusion  in  the  gas 
phase  sre  included  (however,  the  Lewis  number  is  assumed  to  be  unity). 

In  the  original  Hart -Me Cl vie  model  (Ref.  14),  snd  in  the  HFC  model,  the 
burning  rate  is  assumed  to  be  a  function  only  of  the  pressure  and  the 
temperature  gradient  in  the  gas  phase  and  the  latter  is  then  related  to 
the  initial  temperature  of  the  propellant.  By  so  doing,  diffusion  in 
the  gss  phase  is  circumvented.  The  HFC  model  also  takes  into  account 
the  effects  of  a  surface  reaction  between  a  gas  phase  component  and  the 
solid  surface.  However,  the  resultant  response  function  can  easily  be 
reduced  to  what  will  be  referred  to  here  as  the  thermal  case  (essentially 
the  results  derived  by  Bird,  et  al  (Ref.  15)  snd  in  the  same  form  as  the 
Denisor-Baum  results),  or  in  other  words,  the  analysis  without  the  pres¬ 
sure  dependent  surface  reaction.  Throughout  the  remainder  of  this  sec¬ 
tion  the  HFC  model  will  be  considered  as  containing  two  separate  cases, 
one  with  and  one  without  the  surface  reaction. 

If  the  response  functions  of  both  models  (initially  considering  the 
thermal  esse  for  the  HFC  model)  sre  reduced  to  a  common  set  of  nomencla¬ 
ture  and  put  in  the  form  suggested  by  Culick  (see  Appendix  A  of  Ref.  1), 
it  is  found  that  there  is  very  little  difference.  The  DB  response  func¬ 
tion  is 


n 


AB, 


DB 


A3DB  +  I  +  X  '  (A  + 


(4.5) 


and  the  HFC  response  function  is 


1  r 1  fr  A'BHFC  +  n  .,  , 

-  —  j  .  - — - - - - - —  (4.6 

n  p  ^  A'B,  „  +  -7—  +  X  -  (A'  +  1) 

hfc  X 

where  X  is  the  characteristic  root  of  the  energy  equation  for  the  solid 
and  is  equal  to 
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X  =  1/2  [1  +  (1  +  4i)1/2] 


(4.7) 


It  la  apparent  that  the  two  response  functions  are  essentially  the  same 
except  for  a  term  in  the  numerator.  This  term  results  from  the  fact  that 
the  HFC  model  includes  pressure  dependent  surface  reactions.  If  the 
thermal  case  of  the  HFC  model  is  considered  (i.e.,  neglecting  the  sur¬ 
face  reaction)  the  response  function  reduces  to 


r 'Jr  m  _ AEHFC _ 

"  ABHFC+^f +X  "  (A  +  1) 


(4.8) 


Equations  4.5  and  4.8  are  identical,  differing  only  in  the  definition  of 
B,  For  the  two  models  the  definition  of  this  parameter  is 


DB 


* 


and 


(4.9) 


bhfc 


(4.10) 


The  differences  in  the  two  definitions  can  readily  be  traced  to  the 
basic  difference  in  the  two  models  involving  the  form  of  the  burning  rate 
law,  In  the  DB  model  (as  has  been  mentioned)  a  burning  rate  law  is  de¬ 
rived  from  the  energy  and  species  conservation  laws  which  relate  the  burn¬ 
ing  rate  to  the  pressure  and  the  flame  parameters.  In  perturbation  form 
the  term 


n  +  1  +  Ef/2RTf 


is  introduced.  It  has  also  been  pointed  out  that  In  the  HFC  model  the 
burning  rate  is  assumed  to  be  a  function  of  the  pressure  and  the  temper¬ 
ature  gradient  in  the  gas  phase.  The  burning  rate  xs  then  written  as  a 
perfect  differential  in  terms  of  the  appropriate  pressure  and  gradient 
terms.  The  latter  is  then  related  to  the  initial  temperature  of  the 
propellant  which  results  in  the  factor  j/T 
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(where  i  =  ^  ff~)  • 
o 

Thus,  the  only  real  difference  in  the  two  models  being  compared  (consider¬ 
ing  the  HFC  thermal  case  only)  is  in  the  burning  rate  law  that  ia  em¬ 
ployed,  and  for  the  sake  of  numerical  calculations,  the  models  are  identi¬ 
cal.  The  only  difference  lies  in  relating  the  original  parameters  from 
the  respective  models  to  the  value  of  B. 

4.3.2.  Comparing  Calculated  Response  Functions  and  Data 

The  purpose  of  this  section  is  to  compare  the  predictions  of  the 
chosen  instability  models  with  the  results  obtained  from  actual  NAI  data. 
In  order  to  evaluate  Eq.  4.5,  4.6,  and  4.8,  values  for  A  and  B  must  be 
obtained.  Values  of  the  individual  parameters  have  been  measured  or  ob¬ 
tained  from  the  literature  wherever  possible,  and  most  probable  ranges 
of  values  have  been  obtained  where  a  precise  value  is  not  available. 

These  are  tabulated  in  Table  4.1,  with  literature  references  given,  where 
available.  The  corresponding  parameter  values  are  also  tabulated  in 
Table  4.1  where  poasible. 


TABLE  4.1.  Values  of  Parametera  for 
Numerical  Calculation 


Representative 

values 

Values  used  in 
calculations 

Corresponding 

parameter 

values 

T 

0 

300°K 

300°K 

.... 

n 

0. 46a  and  0.49b 

0.50 

.... 

T 

a 

750  to  1,000°K23’24 

900°K 

.... 

E  1 

a 

2225  and  3026  Real 

22,  30  Real 

A  -  8,  11 

j 

0.6  to  0.727’28 

0.65 

bhfc  ■  °-77 

T 

t 

2160d  and  2860°Ke 

2160,  2860°K 

See  Fig.  4.2 

Ef 

10  to  50  Real 

10,  30,  50  Real 

See  Fig.  4.2 

0 

Measured  for  A-35  propellant. 
b  Measured  for  UTAK-TF  propellant  (See  Ref.  22). 
c  Superscripts  refer  to  references. 

d  Calculated  adiabatic  flame  temperature  for  A-35  propellant. 
e  Calculated  adiabatic  flame  temperature  for  UTAH-TF  propellant. 


30 


NWC  TP  4478 


Figures  4.2,  4.3,  and  4.4  show  how  the  nondimen sional  parameters  A 
and  B  vary  as  functions  of  the  physical  parameters.  In  determining  Bpg 
it  was  assumed  that  the  specific  heat  of  the  gas  phase  is  equal  to  that 
of  the  solid,  and  that  the  mean  flame  temperature,  the  mean  surface 
temperature  and  the  burning  rate  exponent,  n,  are  all  constant.  Substi¬ 
tuting  the  values  of  Table  4.1  into  Eq.  4.9,  the  range  of  Bj)g  is  from  0.55 
to  2.0  for  the  two  propellants.  This  range  of  values  reflects  the  un¬ 
certainty  in  the  value  of  the  activation  energy  of  the  flame  which  ia  the 
parameter  having  the  greatest  uncertainty.  This  is  graphically  demon¬ 
strated  in  Fig.  4.2. 


FIG.  4.2.  The  Denison  and  Baum  Gas  Phase  Parameter  B  as  a 
Function  of  E^  With  the  Flame  Temperature  as  a  Parameter. 

From  Fig.  4.3,  it  is  seen  that  the  parameter  &hfc  does  not  vary 
greatly  but  is  close  to  0.8  for  reasonable  values  of  the  parameters. 
The  values  of  the  parameters  in  Table  4.1  result  in  B^fc  “  which 

was  subsequently  used  in  all  the  calculations  for  that  model.  In  Fig 
4.4  A  would  appear  to  have  a  value  on  the  order  of  10.  Combining  the 
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I 


FIG.  4.3.  The  Hart,  Farrell,  Cantrell  Parameter  B  as 
a  Function  of  T  With  the  Temperature  Sensitivity  j 
aa  a  Parameter. 


values  of  the  various  parameters  given  in  Table  4.1  results  in  valuea  of 
approximately  8  and  11  for  the  nonHImensional  parameter  A,  and  these 
values  (for  a  constant  aurface  temperature  of  900°K)  have  been  used  in 
the  calculations  for  both  models. 

Using  the  valuea  of  the  various  parameters  aa  noted  above,  the  real 
part  and  the  phase  of  the  response  function  have  been  calculated  from 
Eq.  4.5,  4.6,  and  4.8  and  are  compared  with  the  data  for  one  of  the  par¬ 
ticular  propellants  in  each  of  the  following  figures.  Thus,  Fig.  4.5 
and  4.6  contain  the  data  of  A-35  propellant  (Ref.  1)  and  the  curvea  cal¬ 
culated  from  the  models  using  the  parameters  that  aie  indicated  in  the 
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FIG.  4.4.  The  Solid  Phase  Parameter  A  Plotted  Versus  Es  With 
the  Surface  Temperature  as  a  Parameter. 


figures.  It  shouid  be  noted  here  that  even  though  the  scatter  in  the 
data  is  considerable,  the  fact  that  a  comparison  can  actuaiiy  be  made 
between  data  of  this  nature  and  the  theoretical  models  is  considered 
significant.  In  comparing  the  two  different  cases  for  the  HFC  modei,  it 
is  apparent  that  the  addition  of  a  pressure-dependent  surface  reaction 
greatly  enhances  the  predicted  response  of  the  propellant. 

Of  the  calculations  for  the  DB  model,  only  the  curve  corresponding 
to  the  fiame-activation  energy  of  30  Kcal/mole  is  included  in  the  two 
piots.  The  curves  for  Ef  •  10  Kcal/moie  resulted  in  responses  that  were 
entirely  too  iow  to  predict  instability,  whiie  the  curves  for  Ej  «  50 
predicted  exponential  instability  (nonsinusoidal  instability;  see  the 
discussion  in  Ref.  16  on  this  subject).  Therefore,  these  curves  were 
not  included  in  Fig.  4.5  and  4.6.  Although  each  of  the  curves  is  quite 
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A- 55  DATA 
E,‘  50  Kc« 

08  Ef  '  50 ,  B ! 0  7 
MFC  (THERMAL  | 
MFC  (PRESSURE) 


FIG.  4.5.  Reaponae  Functiona  aa  a  Function  of  the  Nondimenaional  Fre¬ 
quency.  Si,  and  for  the  Given  Paraoetera  Compared  With  A-35  Data.  HFC 
(thermal)  refera  to  the  HFC  model,  thermal  caae.  HFC  (preasure)  refers 
to  the  HFC  model  including  the  presaure  dependent  aurface  reaction  term 


FIG.  4.6.  The  Response  Functions  for  Es 
Compared  With  A-35  Data. 
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different  from  the  others  it  would  be  impossible  to  judge  from  the  data 
which  set  of  parameters  more  nearly  predicts  actual  behavior.  However, 
the  phase  has  also  been  calculated  for  these  parameters  and  is  compared 
with  the  L*  burner  data  in  Fig.  4.7  and  4.8.  Again,  the  scatter  is  con¬ 
siderable,  but  a  much  more  definite  trend  is  apparent.  On  the  basis  of 
these  latter  figures,  it  can  be  concluded  that  for  both  models  the 
curves  corresponding  to  Es  =  22  Kcal/mole  and,  in  particular,  the  purely 
thermal  part  of  the  HFC  model,  provide  the  better  fit  with  the  data. 

The  results  from  these  latter  two  curves  will  be  used  later  in  a  com¬ 
parison  with  basic  NAI  data. 

! 

Figures  4.9  through  4.12  contain  the  Utah  TF  propellant  data  and 
are  compared  with  the  theoretical  curves.  Two  curves  for  the  DB  model 
have  been  included, ■ corresponding  to  Ef  ■»  30  and  50.  In  Fig.  4.9  and 
4.10  the  value  of  50  provides  the  better  fit  with  the  data.  The  curves 
for  Ej  »  30  are  both  much  lower  than  the  data.  Based  on  these  two 
figures,  it  would  appear  that  the  surface-activation  energy  of  50  in 
both  models  provides  the  better  fit  of  the  data.  Turning  again  to  the 
calculated  values  of  the  phase,  Fig.  4,11  and  4.12  provide  additional 
information  upon  which  a  judgement  msy  be  based.  In  Fig.  4.11,  the 
curve  that  most  closely  fits  the  data  corresponds  to  the  DB  model  with 
Ef  m  30,  but  this  set  of  parameters  has  already  been  rejected  on  the 
basis  of  Fig.  4.8.  From  Fig.  4.12  it  is  again  concluded  that  the  best 
fit  of  the  data  is  obtained  for  Eg  »  22  for  both  models,  utilizing  the 
thermal  part  of  the  HFC  model  (B  »  0.77)  and  Ef  -  50  for  the  DB  model 
(B  -  0.8). 

In  summarizing  the  results  of  this  section,  it  seems  appropriate  to 
choose  the  sets  of  parameters  for  the  calculations  that  result  in  the 
best  fit  of  the  data  for  the  two  models  to  be  coupled  with  NAI  combustor 
theory.  The  only  parameter  that  was  varied  in  the  HFC  model  was  the 
surface-activation  energy  and  the  results  indicate  that  the  value  of 
22  Kcal/mole  (A  »  8)  gives  the  best  results  for  both  propellants.  It 
also  appears  that  the  thermal  case  of  the  HFC  model  is  more  applicable 
than  the  pressure-dependent  case.  This  observation  is  not  apparent  from 
an  examination  of  the  response  function  results  alone  ss  the  real  part 
of  the  two  cases  does  not  vary  much  at  low  frequencies.  However,  the 
phase  is  different  over  the  range  of  frequencies  and  upon  this  basis, 
it  is  possible  to  draw  the  conclusion  that  the  thermal  case  is  more 
applicable.  With  the  DB  model  the  surface-activation  energy  of  22  Kcal/ 
mole  also  gave  the  best  results  for  both  propellants.  This  would  be 
anticipated  because  the  surface-sctivation  energy  is  contained  in  the 
parameter  A  in  exactly  the  same  manner  for  both  models.  However,  for 
the  TP  propellant  a  flame-activation  energy  of  50  Kcal/mole  gave  the  best 
fit  while  the  value  of  30  Kcal/mole  was  best  for  A-35.  A  mechanistic 
argument  explaining  the  difference  between  the  two  propellants  is  beyond 
the  scope  of  this  report  and  will  not  be  attempted  at  this  time.  This 
difference  in  flame-activation  energy  is  not  so  dramatic  when  the  pro¬ 
nounced  differences  of  the  propellant  compositions  are  taken  into 
account. 
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FIG.  4.7.  Theoretical  Phase  Versus  Nondimensional  Frequency  for 
E  -  30,  (A  -  11),  Compared  to  A-35  Data. 
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FIG.  4.11.  rhase  Versus  0  for  Es  =  30,  (A  =»  11),  Compared  to  Utah- 
TF  Data. 
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It  should  be  noted  here  that  the  value  of  B  in  these  two  cases  is 


0.8  and  0.7,  respectively,  while 


5HFC 


was  0.77.  Due  to  the  numerical 


identity  of  the  models  there  is  essentially  one  set  of  nondimensional 
parameters  that  Rive  the  best  fit  for  both  propellants  and  for  both 
models,  i.e.,  A  ■  8,  B  =  0.7-0. 8. 


4.4.  NAI  CHARACTERISTICS  FROM  THEORETICAL  MODELS 

If  the  various  parameters  discussed  in  the  previous  section  were 
known  a  priori,  the  response  function  could  be  calculated  directly  from 
the  models.  This  information  could  then  be  combined  with  Eq.  4.2  to 
give  the  nonacoustic  instability  characteristics  of  a  rocket  motor  or  a 
laboratory  burner. 

In  the  following,  the  h'AI  characteristics  for  the  two  propellants 
studied  will  be  calculated  utilizing  the  models  and  the  apparently  best 
parameters  as  determined  from  Fig.  4.5  through  4.12.  This  will  permit 
the  direct  comparison  of  the  one-dimensional  models  with  experimental 
data  from  an  L*  burner.  Because  of  the  numerical  similarity  between  the 
two  models  only  one  set  of  parameters  was  used  in  the  calculations  for 
each  propellant  (for  A-35:  A  »  8,  B  »  0.77  and  for  Utah  TF:  A  -  8, 

B  -  0.8).  From  Eq.  4.4 


Re  (f^f)  ■  1  +  *Tch 


(4.11) 


and 


tan  f 


ch 


1  +  at 


ch 


(4.12) 


where  if  represents  the  phase  of  the  combustion  relative  to  pressure. 
Combining  these  equations  and  nondimensionalizing  them  with  the  thermal 
wave  time  constant  after  solving  for  t  ^  yields 


ch 


tw 


Re(~7/§)  tan  ♦  Re  (777=)  tan  f 
P  /P  P  IP 


ur 


(4.13) 


tw 


The  right-hand  de  of  the  equation  can  be  completely  determined  theo¬ 
retically  from  tne  models,  thus  relating  the  nondimensional ized  L* 

(t^  *>  L*)  to  the  nondimensional  frequency.  Similar  relationships  have 
been  derived  previously  (Ref.  22),  but  without  considering  the  details 
of  the  gas-phsse  combustion.  Utilizing  the  parameters  for  the  given 
propellants  as  discussed  above,  values  for  the  response  function  and 
phase  as  a  function  of  the  nondimensional  frequency  can  be  substituted 


into  Eq.  4.13  to  give  the  nondimensional ,  L*-frequency  relationship. 

In  Fig.  4.13  theoretical  curves  of  this  nature  are  compared  with  the 
actual  data  for  the  two  propellants.  The  theoretical  curves  are  much 
more  curved  than  one  would  expect  from  the  trend  of  the  actual  data, 
but  other  than  this,  the  agreement  is  quite  reasonable. 

Because  of  the  intimate  relationship  between  the  growth  constant 
and  the  stability  of  the  system,  it  would  be  worthwhile  to  have  a  theo¬ 
retical  relationship  describing  the  growth  constant.  Thia  can  be  de¬ 
termined  readily  by  solving  for  a  from  Eq.  4.12  obtaining 


to 

a  ■  - - 7 

tan  f 


1 


(4.14) 


Evaluating  tc}j  from  Eq.  4.13  and  dividing  through  by  the  frequency 
yields 


a_ 

to 


Re(i^#  - 1 

r  *  fx 

Re  (  , tan  if 
P  /  P 


(4.15) 


Again,  all  of  the  terms  on  the  right-hand  side  of  the  equation  can  be 
determined  from  the  theoretical  models.  Using  the  parameters  noted 
above,  resulta  obtained  from  Eq.  4.15  are  compared  with  actual  data  in 
Fig.  4.14.  Within  the  limits  of  the  data  scatter,  the  theoretical  pre¬ 
diction  is  excellent. 


If  the  growth  constant  is  now  set  equal  to  zero,  the  atability 
limit  for  NAI  is  obtained.  From  Eq.  4.4  and  4.15  this  condition  ia 
aeen  to  be 


1] 


stability  limit 


(4.16) 


From  Fig.  4.14  and  the  equations  for 
ent  that  Eq.  4.16  will  hold  true  for 
quency  (i.e.,  ftcr)  for  any  given  set 
results  into  Eq.  4.13  yields 


the  response  function  it  is  appar- 
a  particular  value  of  the  fre- 
of  parameters.  Substituting  these 


ch 


(taD_i) 


tw 


(4.17) 


Thus,  it  should  be  posalble  to  represent  the  stability  limit  aa  a 
straight  line  paaaing  through  the  origin  on  a  plot  of  tc^  veraus  t1w. 


FIG.  4.13.  Comparison  of  Theory  (Using  the  8est  Parameters  as  Determined  by 
the  Response  Function  and  Phase  Curves)  With  Experiment.  Nnndlmenslonal  L* 
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Figure  4.15  contains  such  a  plot-r-again  comparing  the  results  of  the 
models  with  the  actual  dsta.  The  solid  circles  represent  firings  with 
a  <  0  (i.e.,  in  a  stable  regime  where  developed  oscillations  decayed). 
Considering  that  data  can  only  be  obtained  in  the  unstable  region  and 
up  to  the  stability  limit,  a  line  drswn  along  the  left-hand  edge  of  the 
data  should  represent  the  stability  limit.  For  the  Utah-TF  propellant, 
the  theoretical  results  sgree  quite  well  with  the  data.  However,  for 
the  case  of  the  A-35  data,  it  would  be  virtually  impossible  to  draw  such 
a  line  and  have  it  pass  through  the  origin.  Apparently  something  occurs 
with  this  propellant  which  is  not  accounted  for  by  the  models. 


4.5.  THE  EFFECT  OF  OXIDIZER  PARTICLE  SIZE  ON  NA1 

In  order  to  make  the  models  mathematically  tractable  certain  as¬ 
sumptions  such  as  s  homogeneous  solid  burning  in  one  dimension  with  a 
quasi-steady  gas  phsae  were  made.  In  an  effort  to  assess  the  validity 
of  the  assumption  of  one-dimensional  burning  of  a  homogeneous  solid,  a 
progrsm  was  initiated  to  study  the  effects  of  oxidizer  particle  alze 
distribution  snd  binder  type  on  the  nonacoustic  combustion  instability. 
The  following  sections  present  the  effect  oxidizer  particle  size  has  on 
NA1 . 

4.5.1.  The  Laver  Frequency  Concept 

An  idealized  propel  last  can  be  visualized  sc  consisting  of  oxiui^ei 
particles  stscked  one  on  top  of  another,  surrounded  by  binder.  The 
characteristic  time  to  burn  through  one  particle  is  the  particle  diameter 
divided  by  the  burning  rate  (i.e.,  T  •  D/r) .  If  the  particles  are  in¬ 
deed  stacked  in  layers  then  this  characteristic  time  will  be  repested  at 
a  frequency  of  one  cycle  per  particle,  or  f  -  1/t  •  r/D  (Fig.  4.16).  If 
the  entire  surface  of  a  burning  propellant  regressed  so  that  each  oxi¬ 
dizer  particle  was  exposed  at  the  ssme  time  and  burned  "in  phase"  with 
enough  other  oxidizer  particles,  then  this  layer  frequency  would  become 
a  measurable  phenomenon,  contributing  to  combustion  instability.  This 
phase  correlstion3  concept  is  schematically  pictured  in  Fig.  4.17. 

In  sn  actual  propellsnt  the  oxidizer  psrticles  are  rsndomly  dis¬ 
tributed  and  therefore  the  phase  correlation  of  the  entire  surface  would 
be  somewhat  difficult  to  schieve  and  would  most  likely  result  in  8  fre¬ 
quency  less  than  noted  above.  That  is  to  say, 

f  -  K/t  -  Kr/D  (4.18) 


3  This  "phase  correlation"  concept  was  first  proposed  by  Price  (Ref.  29) 
in  connection  with  the  shedding  of  sgglomerated  metal  from  a  propellant 
surface. 
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FIG.  4.15.  Comparison  of  Theories  With  Experiment.  Stability  limit-chamber 
residence  time  versus  thermal  wave  relsxation  time.  To  the  right  of  the 
lines  is  the  unstable  domain. 
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FIG.  4.16.  The  Lsyer-Frequency  Concept — The  Frequency  at  Which 
a  Combustion  Wave  Burns  Through  Layered  Particles  Is  Directly 
Proportional  to  the  Burning  Rate  and  Inversely  Proportional  to 
the  Diameter. 


Phase  correlated 
surface  with 
particles  burning 
in  phase 


FIG.  4.17.  Phase  Correlation — If  the  Burning  Surface  Arranges 
Itself  Such  That  Sufficient  Oxidizer  Particles  are  Exposed 
Simultaneously,  Then  the  Surface  Is  Said  to  be  "Phase  Corre¬ 
lated",  and  Measurable  Oscillations  will  Occur  at  the  "Layer 
Frequency". 
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where  <  1  due  to  the  fact  that  the  oxidizer  particlea  are  not  atacked 
uniform../  one  on  top  of  the  other  but  are  randomly  oriented.  The  poasi- 
bility  of  having  a  phaae  correlated  aurface  can  be  explored  by  examining 
the  surface  of  a  quenched  propellant.  Figure  4.18  is  a  photomicrograph 
of  the  surface  of  a  propellant  that  extinguished  Itself  while  chuffing. 
The  protruding  surfaces  of  the  oxidizer  particlea  are  readily  apparent 
and  It  is  not  unrealistic  to  imagine  that  the  aurface  was  indeed  phase 
correlated  to  aome  degree.  The  extent  to  which  a  surface  must  be  "corre¬ 
lated"  in  order  for  instability  to  be  obaervable  la  a  queation  that  could 
be  pursued  in  the  future. 


FIG.  4.18.  Photomicrograph  of  an  Extinguished  Sample  of  A-155 
Propellant.  The  large,  protruding  crystala  had  an  initial, 
mean  diameter  of  400p. 


It  was  pointed  out  in  Section  4.1  that  a  dual  frequency  nature  had 
been  exhibited  by  two  different  propellanta,  and  that  thia  duality  had 
been  attributed  to  the  oxidizer  particle  aize.  According  to  the  layer 
frequency  concept  outlined  above,  any  propellant  containing  a  bimodal 
blend  of  oxidizer  should,  under  the  right  conditions,  display  this  dual¬ 
frequency  behavior.  To  check  this  postulate  a  series  of  bimodal  AP  pro¬ 
pellants  was  made  having  252  polyurethane  binder  and  752  solids  (12 
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carbon  black  was  used  to  reduce  the  translucence  of  the  propellant) . 

Part  of  the  AP  in  the  propellants  was  a  coarse  oxidizer  and  part  was 
finer  oxidizer,  with  the  oxidizer  particle  size  distribution  for  each 
propellant  given  in  Table  A. 2.  The  AP  that  was  used  was  the  special 
spherical  AP  produced  by  the  American  Potash  and  Chemical  Corporation, 
and  all  blends  except  the  15p  were  screened  to  give  either  +18  or  +25% 
of  the  mean  particle  diameter.  Chemically,  the  propellants  were  the 
same;  only  the  physical  nature  was  changed  by  varying  the  oxidizer  parti¬ 
cle  sizes.  In  the  remainder  of  this  section  the  results  obtained  with 
these  propellants  and  their  relationship  with  the  layer  frequency  will  be 
discussed. 


TABLE  A. 2.  Bimodal  Propellant  Compositions. 


Ingredients  and  weight  percent 


Propellant 

designation 

1 

Ammonium  Perchlorate 

Binder 

Other 

A-1A6 

i 

■37.5% 

15(i 

37.5% 

80p 

25%  polyurethane 

A-1A8 

37.0% 

15(1 

37.0% 

200(1 

25% 

do 

1% 

carbon  black 

A-1A9 

37.0% 

90y 

37.0% 

600y 

25% 

do 

1% 

do 

A-151 

37.0% 

A5p 

37.0% 

200p 

25% 

do 

1% 

do 

A-155 

37.0% 

A5p 

37.0% 

AOOu 

25% 

do 

1% 

do 

A-156 

51.8% 

15u 

22.2% 

200y 

25% 

do 

1% 

do 

A-157 

22.2% 

15u 

51.8% 

200u 

25% 

do 

1% 

do 

A-158 

37.2% 

90y 

37.5% 

600(i 

25% 

do 

A-159 

36.0% 

15(i 

36.0% 

200u 

25% 

do 

1% 

carbon  black 

1% 

n-butyl  ferro- 

cene 

A-160 

36.0% 

A  5u 

36.0% 

A00M 

25% 

do 

1% 

carbon  black 

1% 

n-butyl  ferro- 

cene 

A. 5. 2.  The  Non-one-dlmenslonal  Nature  of  Propellant 

The  models  available  (Section  A. 3)  to  predict  the  general  behavior 
of  a  given  propellant  all  assume  that  the  propellant  is  a  homogeneous 
solid  burning  in  a  one-dimensional  manner.  The  degree  of  validity  of 
these  assumptions  may  be  investigated  by  examining  the  heterogeneous, 
non-one-dimensional  nature  of  the  propellant  with  respect  to  thermal 
wave  penetration  and  oxidizer  particle  size.  Figure  A. 19  is  a  schematic 
showing  the  relative  dimensions  of  typical  oxidizer  particles  and 
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FIG.  4.19.  Thermal  Waves  In  Heterogeneous  Propellants— The 
Circles  Represent  AF  Crystals  (Drawn  to  Scale)  and  the  Dashed 
Lines  Represent  the  Temperature  Profile  and  the  Distance  thst 
the  Profile  Penetrates  Into  the  Crystal.  The  upper  profile 
was  calculated  from  the  properties  of  AP  alone  end  the  bottom 
one  Is  based  on  the  average  propellant  properties. 


temperature  profiles  for  two  different  burning  rates.  The  thermal  pro¬ 
files  which  have  been  Included  are  based  on  the  thermal  properties  of 
AP  alone  (top  profile)  and  on  the  average  properties  of  the  propellant 
(bottom  profile).  At  low  pressure  and  low  burning  rate,  the  thermal 
wave  thickness  Is  quite  large,  and  for  small  oxidizer  particle  elzes 
the  propellant  can  be  considered  hosogeneous  and  burning  may  essentially 
be  considered  one-dimensional.  However,  at  a  burning  rate  of  0.1  ln/sec, 
the  thermal  wave  thickness  Is  the  same  order  of  magnitude  as  the  oxidizer 
particle  sizes  used  In  this  atudy.  Hence,  the  homogeneoue  and  one- 
dlmenslonal  assumptions  necessary  to  the  analytical  modele  would  appear 
to  be  tenuous  under  these  conditions,  snd  et  the  same  time  these  ere  the 
conditions  where  the  layer  frequency  would  most  likely  be  observed.  If 
the  thermal  wave  thickness  just  matches  the  particle  size,  then  as  the 
particle  burn3  the  combustion  wave  will  soon  be  passing  into  preheated 
AP  and  the  local  burning  rate  of  the  particle  will  increase  as  the 
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particle  burns  out.  This  rspld  burning  of  the  remains  of  an  oxidizer 
particle  has  been  observed  in  movies  such  as  those  described  in  Section 
3.  It  can  readily  be  imagined  that  if  such  a  phenomenon  could  be  "phase 
correlated"  over  the  entire  propellsnt  surface,  pressure  oscillations 
would  be  produced. 

4.5.3.  Results  From  Blmodal  Propellants 

From  the  above  discussions  the  greatest  departure  from  one- 
dimensional  theory  would  be  expected  st  high  values  of  burning  rate 
(smaller  thermal  wave)  and  for  large  oxidizer  psrtlcle  sizes.  This  Is 
indeed  the  case  as  illustrated  by  the  trends  of  the  data  plotted  in 
Fig.  4.20.  Test  data  as  well  as  the  burning  rates  of  the  various  pro¬ 
pellants  are  tabulated  in  Appendix  A.  The  cross-hatched  parabolic  band 
represents  the  prediction  of  one-dimensional  theory,  and  the  data  are 
for  the  various  propellants.  At  low  vslues  of  the  burning  rate, 
r  <  .09  in/sec,  the  dsta  agree  moderately  well  but  at  higher  burning 
rates  (r  >  .10  in/sec),  the  experimental  values  deviate  quite  signifi¬ 
cantly  from  the  one-dimensionsl  theory. 

Whereas  Fig.  4.20  is  a  composite  plot,  Fig.  4.21-4.26  present  the 
data  for  each  individual  propellsnt.  In  the  figures  the  shaded  areas 
bounded  by  straight  lines  represent  the  mapping  of  the  equation  of 
f  »  r/P  for  the  particle  size  distribution  used  in  the  formulation  of 
that  propellant.  These  shaded  areas  show  where  the  layer  frequency 
should  show  the  strongest  effect  (for  K  -  1). 

Figure  4.21  prosents  the  dats  for  the  unimodal,  90u  AP  propellant, 
A-35.  For  this  propellant  the  data  agree  well  with  both  the  layer 
frequency  concept  (the  line  labeled  85p)  and  the  one-dimensional  theory 
as  the  two  are  coincident  for  the  range  of  burning  rates  available  for 
testing. 

Figure  4.22  presents  the  data  for  the  blmodal  (90  snd  600u)  pro¬ 
pellant  A-149.  Since  the  large  particles  of  AP  tend  to  decrease  the 
burning  rate,  it  is  to  be  expected  that  the  layer  frequency  concept  (for 
the  90p  AP)  and  the  one-dimensional  theory  are  mutually  applicable  for 
this  propellant.  Apparently  the  layer  frequency  corresponding  to  the 
600u  AP  was  so  low  that  it  was  unobservahle. 

The  one-dimensional  theory  parabola  deviates  quite  significantly 
from  the  experimental  results  for  the  A-151  propellant  (45  snd  200u)  as 
shown  in  Fig.  4.23.  As  was  mentioned  earlier  for  a  burning  rate  of 
0.10  in/sec,  the  thermal  wave  thickness  for  the  propellsnt  is  approxi¬ 
mately  45u  which  coincides  with  the  fine  oxidizer  particle  size  for  this 
propellant.  Hence,  non-one-dimens ional  behavior  was  expected  for  this 
range  of  conditions.  It  should  be  noted  here  as  well  as  with  the  other 
propellants  that  the  range  of  conditions  for  which  data  may  be  obtsined 
is  very  limited  in  an  L*  burner.  At  low  burning  rates  where  the  steady- 
state  thermal  wave  penetrates  deeply  into  the  solid  there  are  ignition 
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FIG.  4.20.  Frequency  Veraue  Burning  Rate  for  the  Propellanta 
Teeted.  The  data  for  the  reference,  unimode  1  propellant  are 
Indicated  by  the  ahaded  area  and  the  predicted  behavior  (baeed 
on  one-dlnenalonal  nodela)  le  Included. 
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FIG.  4.21.  Frequency  Versus  Burning  Rate  for  the  Reference 
Propellant,  A-35.  The  straight  line  represents  the  layer 
frequency  concept  based  on  the  median  of  the  particle  dis¬ 
tribution  of  the  AP. 


problems  In  establishing  this  profile  and  Igniting  the  propellant.  A  so, 
there  are  almllar  problems  as  the  low  pressure  deflagration  limit  Is 
reached.  At  high  burning  rates,  data  acquisition  Is  restricted  by  the 
stability  limit  of  the  L*  burner,  and  by  the  physical  limit  of  the  burner 
(i.e.,  the  smaller  values  of  L*  that  would  be  below  the  stability  limit 
are  physically  Impossible  to  attain  In  the  burner). 

Figures  4.24  and  4.25  are  quite  similar  as  might  be  expected,  due 
to  the  similarity  in  the  propellant  compositions.  Both  contain  15u  AP 
and  both  oscillated  at  frequencies  two  to  three  times  what  normally 
would  have  been  expected.  It  appears  as  though  there  Is  a  transition 
occurring  from  one-dlmenslonal  behavior  at  low  burning  rate  to  a  layer 
frequency  behavior  at  the  higher  burning  rate. 

In  Fig.  4.26,  propellant  A-155  not  only  exhibits  the  higher  trend 
in  frequencies  but  also,  under  certain  conditions,  oscillated  at  very 
low  frequencies.  The  low  frequency  corresponds  very  closely  to  the 
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F1C.  4.22.  Frequency  Veraus  Burning  Rate  for  A-149. 


layer  frequency  cf  the  larger  oxidizer  partlclee.  This  represents  the 
dual  frequency  nature  thet  wae  aought  at  the  outaet  of  tha  atudy. 
Apparently  e  sufficient  number  of  the  large  oxidizer  pertlclea  ware 
burning  in  pheae  with  eech  other  to  cauae  en  oscillation  of  sufficient 
magnitude  that  it  could  be  meaeured.  Figure  4.27  la  a  preesure-tlae 
trece  of  one  of  the  rune  where  the  low  frequency  oadlletlon  occurred. 

The  chuffing  end  oscllletlons  thet  occur  at  the  initial  part  of  the  test 
ere  at  a  frequency  of  52  cpe  while  the  low  amplitude  oaclllatlona  later 
in  the  teat  ere  at  e  frequency  of  2.7  cps.  Thle  la  the  moet  concluelve 
evidence  avelleble  thet  the  layer  frequency  concept  la  velld  end  producea 
en  oheervable  affect  under  certeln  conditiona.  Subsequent  test  records 
have  elao  exhibited  thle  dual  frequency  phenomenon. 
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FIG.  4.23.  Frequency  Versus  Burning  Rate  for  A-151. 

4.6.  THE  RELATIONSHIP  BETWEEN  THEORETICAL  MODELS 
AND  THE  LAYER  FREQUENCY  CONCEPT 

In  order  to  put  the  results  of  this  study  In  perspective  with 
Section  4.3,  the  data  have  also  been  reduced  to  give  the  response  func¬ 
tion  and  the  nondlmenslonal  frequency.  Figure  4.28  contains  the  re¬ 
sults  of  the  data  for  four  of  the  blmodal  propellants  and  compares  these 
with  results  obtained  using  the  Denison  and  Baum  model.  The  parameters 
used  In  the  theoretical  calculation  were  determined  for  A- 35  propellant 
In  the  manner  outlined  In  Section  4.3  and  also  by  the  method  recently 
proposed  by  Beckstead  and  Culick  (Ref.  30).  Utilizing  this  latter 
technique  to  determine  the  parametet*  necessary  to  fit  the  data  for  the 
blmodal  propellants  results  In  parameters  that  are  totally  Impractical 
(l.e. ,  A  <  0).  This  again  Indicates  that  these  propellants  oscillate 
In  a  manner  that  cannot  be  predicted  from  one-dlmenalonal  theories. 

Another  way  of  looking  at  the  data  la  In  e  stability  plot  similar 
to  Fig.  4.15.  Utilizing  the  data  for  A-155,  which  seemed  to  show  the 
most  dramatic  dual  frequency  results.  Fig.  4.29  waa  constructed.  The 
data  for  o  >  0  should  lie  to  the  right  of  the  stability  limit  line,  the 
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FIG.  4.24.  Frequency  Vereus  Burning  Rate  for  A-146. 
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Frequency  Versu.  Burning  Rate  for  A-14B 


FIG.  4.28.  Response  Function  Vereue  Nondlnensional  Frequency. 
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FIG.  4.29.  Stability  Limit  Plot  of  tch  Versus  Ttw  for  A-155. 
Dst.a  t0  the  right  of  the  lines  are  In  the  unstable  region. 


a  -  0  points  should  lie  on  the  line,  snd  the  a  <  0  lines  should  be  to  the 
left  of  the  limit  In  the  stsble  domain.  In  the  figure  all  of  the  points 
with  a  >  0  lie  to  the  right  of  the  solid  line  that  Intersects  the  ab¬ 
scissa  et  itw  -  20  milliseconds.  The  majority  of  the  date  lying  to  the 
left  of  this  llae  are  the  date  et  very  low  frequency.  This  seems  to 
lndlcete  the  possibility  of  two  different  stability  limits,  one  corre¬ 
sponding  to  the  very  low  frequency  date  end  the  second  for  the  remainder 
of  the  date.  It  will  be  recalled  from  Section  4.4  that  the  theories 
predict  e  atrelght  line  through  the  origin  es  the  eteblllty  limit.  In 
Fig.  4.29  the  very  low  frequency  dete  seem  to  give  a  stsblllty  limit 
that  would  egree  with  theories  and  yet  the  frequencies  ef  the  oscil- 
letlons  ere  much  lower  than  would  be  predicted.  This  anomaloue  behavior 
nteds  further  lnveetlgetlon  In  order  to  underetend  the  relation  between, 
and  the  reletlve  effect  of,  the  leyer  frequency  end  one-dlmensionel 
theories. 
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From  the  ahove  discussion  it  is  apparent  that  the  one-dimensional 
theories  available  at  present  are  not  adequate  to  predict  the  unstable 
behavior  of  propellants  In  a  quantitative  manner.  Nor  is  it  possible 
to  do  so  with  a  simple  mechanistic  argument  such  as  the  layer  frequency 
concept.  However,  by  combining  aspects  of  each  of  these  a  greater  under¬ 
standing  of  the  problem  may  possibly  be  attained.  In  the  following  dis¬ 
cussion  it  should  be  kept  in  mind  that  the  layer  frequency  concept  re¬ 
sults  in  a  range  of  frequencies  under  which  oscillations  occur  but  gives 
no  indication  as  to  the  magnitude  of  the  oscillations,  whereas  the  one¬ 
dimensional  models  give  both  the  frequency  and  an  indication  of  the 
magnitude  of  the  oscillations  (in  the  form  of  the  response  function). 
Therefore,  the  two  concepts  cannot  be  combined  in  an  exact  sense,  but  a 
qualitative  combination  of  the  results  is  very  informative.  A  plot 
similar  to  Fig.  4.20-4.26  has  been  constructed  in  which  the  results  of 
the  one-dimensional  models  can  be  readily  included.  Figure  4.30  is  a 
log-log  plot  of  the  burning  rate  and  the  frequency  with  the  layer  fre¬ 
quency  lines  corresponding  to  A-155  propellant  included  as  well  as  the 
line  representing  the  maximum  in  the  response  function  for  A  “  14  and 
B  ■  0.8  (see  Fig.  4.28).  If  a  third  axis  extending  outward  from  the 
plane  of  the  paper  is  considered  to  represent  the  real  part  of  the  re¬ 
sponse  function,  then  Fig.  4.30  takes  the  form  of  a  contour  map.  On  the 
map  there  will  be  three  ridges,  two  corresponding  to  the  layer  fre¬ 
quencies,  and  the  third  corresponding  to  the  maximum  in  the  one¬ 
dimensional  theory.  By  taking  a  section  of  the  contour  at  a  constant 
burning  rate  (r  ■  0.09  in/sec),  a  modified  form  of  the  response  function 
can  be  postulated  as  in  Fig.  4.31.  The  dashed  lines  represent  the  approx¬ 
imate  effect  that  could  be  caused  by  the  layer  frequency,  and  it  has  been 
assumed  that  the  two  effects  would  be  additive.  The  magnitude  of  these 
ridges  due  to  layer  frequency  effects  are  unknown  before  a  series  of 
tests  are  conducted,  but  Fig.  4.32  may  be  used  to  obtain  the  frequency 
that  might  be  expected,  given  the  oxidizer  particle  size  and  a  particular 
burn  rate.  This  figure  also  shows  how  this  result  would  differ  from  that 
expected  from  a  simple  one-dimensional  analysis,  such  as  the  Denison  and 
Baum  model. 


4.7.  SUMMARY 

Two  prevalent  analytical  models  of  solid  propellant  combustion  in¬ 
stability  which  are  relevant  to  the  low-frequency  regime  have  been  re¬ 
duced  to  a  common  notation  and  found  to  be  numerically  identical.  Non¬ 
acoustic  data  for  several  propellants  have  been  reduced  to  a  format  common 
with  the  models,  thus  allowing  a  quantitative  comparison  between  theoreti¬ 
cal  analyses  and  experimental  data  of  this  nature,  The  possibility  of 
utilizing  combustion  dynamics  models  to  predict  nonacoustic  instability 
behavior  has  been  explored:  only  qualitative  agreement  between  experi¬ 
mental  data  and  model  predictions  has  been  found,  and  even  then  several 
large  discre  ancles  have  been  discovered.  One  such  discrepancy,  the 
effect  of  oxidizer  particle  size,  has  been  explored  in  detail.  A  very 


63 


t  (in/*«c) 


FIG.  4.30.  Burning  Rate  Versus  Frequency  for  A-155  Propellant 
With  Shaded  Areas  Representing  the  Layer  Frequency  Regimes. 

The  heavy  solid  line  is  the  maximum  of  the  response  function 
calculated  from  the  Denison  and  Baum  model. 


simple,  phenomenological  argument  relating  frequency  of  oscillation  to 
the  particle  size  of  the  oxidizer  has  been  postulated  and  experimentally 
examined.  This  was  accomplished  by  formulstlng  s  series  of  propellants 
containing  bimodsl  oxidizer  blends  snd  testing  them  in  the  L*  burner. 

The  results  obtained  can  only  be  explained  by  considering  aspects  con¬ 
tained  in  the  models  as  well  as  aspects  of  the  layer  frequency  concept 
which  is  related  to  the  oxidizer  particle  size. 
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FIG.  4.31.  The  Modified  Response  Function  With  Deviations 
Included  (the  Dashed  Lines)  to  Account  for  Oxidizer  Particle 
Size. 
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FIG.  4.32.  Burning  Rate  Versus  Frequency  With  the  Layer  Fre¬ 
quency  Lines  Included  for  Typical  AP  Sizes  and  the  Theoretical 
Maximum  Calculated  From  a  Denison  and  Bsum  Type  Model. 
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5.  ACOUSTIC  COMBUSTION  INSTABILITY 


Experimental  equipment  has  been  devised  which  permita  the  investi¬ 
gation  of  acoustic  instability  in  the  low  frequency  range  characteristic 
of  very  large  rocket  motors  (Ref.  31,  32,  33  and  34)  snd  recent  modifi¬ 
cations  to  the  equipment  permit  control  of  the  experimental  parameters 
and  acquisition  of  quantitative  information  Including  acoustic  response 
function  data  (Ref.  1).  A  cutawsy  diagram  of  the  burner  used  in  this 
lnvestlgstlon  is  shown  in  Fig.  1. 

Activity  during  the  present  report  period  included  assessment  of 
stability  characteristics  of  several  metallized  and  nonmetalllzed  pro- 
pellanta,  initial  tests  to  determine  the  acoustic  response  function  of  a 
nonmetalllzed  polyurethane  propellant  composition,  tests  with  pressed 
AP-A1  pellets  which  indicate  that  these  spontaneously  unstable  formu¬ 
lations  will  couple  with  and  drive  acoustic  wave  modes,  and  an  investiga¬ 
tion  of  nonlsentroplc  phenomena  in  low  frequency  acoustic  instability  of 
solid  propellants. 


5.1.  TESTS  TO  DETERMINE  RANGE  OF  UNSTABLE  BEHAVIOR 

Tests  were  conducted  to  determine  unstable  combustion  behavior  of 
two  compositions  for  big  booster  applications  (ANB  3105  and  ANB  3254), 
a  metallized  and  a  nonmetalllzed  AP-polyurethsne  propellant  (A-147  and 
A-35),  and  a  nonmetalllzed  AP-PBAN.  composition  (A-13).  Results  of  the 
big  booster  propellsnt  tests  are  dlscuased  in  Appendix  B. 

An  investigation  previously  conducted  with  a  family  of  metallized 
AP-polyurethane  propellants  (Ref.  .'.)  to  determine  the  range  of  unatable 
behavior  was  expanded  to  include  a  new  composition.  Previous  tests  on 
A-139  snd  A-91  propellants  revealed  that  A-139,  although  unstable  over  a 
much  wider  range  of  pressure  and  frequency  than  A-91,  was  much  hsrder  to 
ignite  (information  on  propellant  compositions  appears  in  Table  5.1). 
Since  the  AP  psrtlcle  size  was  the  only  parsmeter  vsried  between  A-139 
and  A-91  it  was  expected  that  a  propellant  with  an  intermediate  oxidizer 
size  might  offer  a  wider  range  of  instability  than  A-91  and  at  the  same 
time  be  easier  to  ignite  than  A-139.  A  50-pound  mix  of  A-147  was  made 
using  45p  spherical  AP.  Tests  to  determine  the  range  of  instability 
were  conducted  and  the  results  are  shown  in  Fig.  5.2  with  A-139  and  A-91 
included  for  comparison.  The  rsnge  of  instability  for  A-147  appeared  to 
be  less  than  that  for  A-91  and  ignltablllty  of  A-147  did  not  appear  to 
be  greatly  improved  over  A-139;  hence,  further  tests  with  A-147  are  not 
planned. 
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TABLE  5.1.  Compositions  Used  in  Low  Frequency 
Acoustic  Instability  Research 


Designation 

of 

composition 

Binder 

Oxidizer 
(spherical  AP) 

Alumint 

(spheric 

m 

al) 

Other 

Type 

wt  X 

Mesn  dls. 

Wt  X 

Mesn  dis. 

wt  x 

Type 

Wt  X 

A-13 

PBAN 

24.0 

80  u 

76.00 

— 

— 

— 

A-35 

Polyur- 

25.0 

80  y 

75.00 

— 

1 

— 

— 

ethane 

A-91 

do 

25.0 

80  y 

67.00 

!  5  R 

8 

— 

— 

A-139 

do 

25.0 

15  v 

67.00 

5  y 

8 

— 

— 

A-147 

do 

25.0 

45  y 

67.00 

5  u 

8 

— 

— 

A-152 

do 

25.0 

80  V 

74.60 

— 

1  * 

NaCl 

0.4 

20.1 

5  y 

23.85 

— 

NaCl 

0.05 

PU 

do 

— 

— 

— 

45  y 

56.00 

■ 

J-30 

— 

1 

80  U 

94.00 

5  y 

Cu0202 

2.00 

J-81 

— 

80  y 

88.00 

15  u 

■ 

CuO 

8.00 

Two  unmetallized  propellants,  A-13  snd  A-35  (Table  5.1),  which  have  been  ex¬ 
tensively  tested  in  the  1.5-inch  T-burner,  were  chosen  for  study  in  the  low  fre¬ 
quency  range  in  the  5.5-inch  burner.  Testing  of  these  two  compositions  to  dste 
has  been  primarily  to  determine  the  range  of  unstable  combustion  behavior. 

Results  are  shown  in  Fig.  5.3. 

A  notable  feature  of  tests  with  A-13  was  the  simultaneous  presence  of 
several  sxlsl  acoustic  modes  during  tests  when  the  fundamental  mode  frequency  was 
below  90  cps.  Validity  of  the  response  function  determination  depends  on  pres¬ 
ence  of  s  single  acoustic  mode  s.  ;  >■  st  A-13  response  function  tests  sppesr  to  be 
limited  to  frequencies  above  sbout  90  cps.  Tests  with  A-35  over  s  range  of  fre¬ 
quencies  similar  to  the  A-13  tests  did  not  reveal  presence  of  modes  higher  than 
the  fundamentsl  axisl  mode  snd  there  does  not  appear  to  be  s  combustion- imposed 
lower  frequency  limit  ss  with  A-13. 
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MAH  POITUKTHAHC 


FIG.  5.3.  Ranga  of  Instability  Behavior 
for  Two  Unmetallized  Propallants. 


5.2.  PRESSED  PELLET  INSTABILITY 

Spontaneous  combustion  oscillations  observed  with  AP-A1  pressed 
pallets  have  been  previously  reported  (Ref.  1).  Investigations  with 
metallized  propellents  which  appeared  to  hava  a  preferred  frequency  be¬ 
havior  hava  also  bean  published  (Ref.  35),  but  the  correlation  of  s  pre¬ 
ferred  frequency  charecteristic  with  metsl  combustion  Jn  propellant  com¬ 
positions  has  been  indirect.  Thus,  preseed  AP-A1  materials  appear  to 
offer  e  means  for  conducting  s  definitive  study  of  the  relationship  be¬ 
tween  metsl  combustion  snd  preferred  frequency  combustion  oscillations. 

The  relevance  of  pressed  pellet  combustion  oscillations  to  acoustic 
instability  has  not  been  tested  until  recently.  At  leset  two  aspects 
are  Involved  which  determine  whether  epontanaous  combustion  oscillations 
will  excite  an  acoustic  wava  mods  oscillation:  (1)  whether  the  combustion 
oscillations  generete  e  pressure  perturbation,  end  (2)  whether  the  com- 
buetlon  oscillations  ere  capable  of  being  coupled  with  flow  disturbances 
in  a  manner  which  will  result  in  amplification  of  sn  acoustic  mode 
(Ref.  36). 

Initial  tasts  to  determine  the  extent  to  which  preesed  AP-A1  ma- 
terlale  exhibit  ecouetlc  insteblllty  have  been  conducted  with  J-30  (see 
Teble  5.1)  in  the  5.5-lnch  burner.  The  probleme  of  ettelnlng  lerge 
pellet  size  and  of  swuntlng  the  pellets  without  the  uee  of  organic 
blndere  or  potting  compounds  which  inhibit  spontaneous  instability, 
were  eolved  by  utilizing  sn  exletlng  die  of  approximately  1-1/2-in. 
diameter  and  erranglng  the  pellets  in  en  errsy  such  that  seven  pellete 
ware  acconnodeted  in  each  end  of  the  burner.  Mounting  wss  accomplished 
by  pressing  the  pellets  into  a  metsl  disc  (Fig.  5.4).  The  metal  disc- 
pellet  aseembly  wse  then  held  in  plsce  in  the  bottom  of  the  burner  end- 
cep  by  s  backing  of  slllcona-gless  adhesive  tape.  Ignition  wss  accom¬ 
plished  by  use  of  the  eame  paint-on  msterlsl  used  for  propellants. 
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FIG.  5.4.  Pressed  Pellet  Mounting 
Arrsngenent  for  the  5.5-inch  Burner. 


Test  results  Indicated  that  coupling  did  occur  and  that  the  pressed 
pellet  material  la  a  good  medium  for  atudy  of  the  role  of  metsl  behavior 
in  producing  preferred  frequency  Instability.  Further  testing,  however, 
will  be  more  resdlly  accomplished  by  choice  of  a  pressed-pellet  formu¬ 
lation  that  exhibits  spontaneous  oscillatory  behavior  st  lover  pressures 
than  J-30.  J-81  (Table  5.1)  appears  to  be  a  suitable  candidate,  a-  noted 

In  Ref.  1. 


5.3.  INVESTIGATION  OF  NON1SENTROPIC  PHENOMENA 
5.3.1.  Introduction 

Knowledge  of  the  response  of  propellant  combustion  processes  to  flow 
perturbations  Is  vital  to  any  quantitative  description  of  acoustic  Insta¬ 
bility.  The  acoustic  admittance  appears  as  a  boundary  condition  In  so¬ 
lution  of  the  wave  equation  and  It  describes  the  flow  of  energy  Into  or 
out  of  the  acoustic  wave.  The  admittance  la  defined  as 


Y  = 


(5.1) 


71 


t 


NWC  TP  4478 


where  v'  is  s  velocity  perturbation  normal  to  the  propellant  surface  and 
p’  is  a  pressure  perturbation.  When  written  in  terms  of  variables  re¬ 
lated  to  propellant  combustion,  the  admittance  becomes  (Ref.  34) 


where  v  is  the  mean  gas  velocity,  'p  la  the  mesn  pressure,  u/c  la  the  ratio 
of  mssa  perturbation  to  pressure  perturbation,  and  o/e  is  the  ratio  of 
density  to  pressure  perturbation.  It  Is  assumed  that  the  above  quanti¬ 
ties  are  calculated  at  a  plane  near  the  propellant  surface  but  at  a  loca¬ 
tion  outside  the  chemical  reaction  zone.  It  Is  further  assumed  that  gsa 
velocity  and  pressure  fluctuations  aro  described  by  s  standing  acoustic 
wave  mode,  that  the  flame  zone  is  of  a  thickness  small  compared  to  other 
dimensions  in  the  system,  snd  that  the  gss-phase  reaction  times  are  short 
compared  to  a  period  of  acouatic  oscillation  snd  to  the  transit  time  of 
sn  element  of  gas  through  the  combustion  zone.  It  should  be  noted  that. 

In  geneisl,  Y  is  a  complex  quantity.  However,  the  real  part  of  the  com¬ 
plex  admittance  determines  whether  or  not  an  acouatic  wsve  wiil  be  ampli¬ 
fied;  the  out-of-phase  (or  Imaginary)  component  of  the  combustion  re¬ 
sponse  cannot  amplify  a  pressure  perturbation. 

The  ratio  u/c  ia  commonly  referred  to  aa  the  "response  function"  aid 
considerable  theoretical  snd  experimental  Investigation  has  been  devoted 
to  determining  its  magnitude.  Values  of  the  response  function,  according 
to  T-burner  results,  typically  fall  in  the  range  between  unity  snd  S. 

Less  emphasis  hsa  been  put  on  evaluating  o/e  since  theoretlcsl  studies 
indicate  that  It  Is  of  the  order  of  unity  in  magnitude.  It  Is  this  ratio 
which  is  sensitive  to  the  nature  of  the  gas  processes  in  the  combustion 
zone.  If  the  combustion  zone  behaves  lsentropically,  o/c  -*■  while  it 
approaches  unity  If  the  combustion  zone  behaves  isothermally.  This  latter 
situation  is  often  referred  to  as  the  "zero-frequency  caae".  Actual 
values  of  this  rstlo  probably  lie  between  these  two  extremes,  varying 
according  to  conditions. 

In  spite  of  the  reduced  importance  of  o/c  to  determination  of  the 
admittance,  Interest  in  experlmentel  investigation  of  conditions  under 
which  ncnisentropic  behavior  sppeera  has  persisted.  Attempts  to  verify 
nonlsentioplc  behavior  during  acoustic  instability  have,  until  the  present, 
been  negative  or  Inconclusive.  In  theory  the  Isothermal  case  assumes 
thst  the  flame  temperature  In  an  oscillating  environment  Is  constant. 

Under  typical  rocket  motor  operating  conditions  the  flame  temperature  is 
relatively  insensitive  to  pressure  t>v  the  above  assumption  would  appear 
to  be  reasonably  valid.  Under  these  conditions,  successive  elements  of 
gae  enter  the  ecoustlc  field  with  different  values  of  entropy  which  Is 
given  by  the  equation  (Ref.  37) 


a  -  a  +  cR(sln  {wt-x]), 
o 


(5.3) 
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where  a0  is  a  reference  value  of  entropy,  e  is  the  ratio  of  maximum  pree- 
aure  fluctuation  to  mean  preaaure,  R  ia  the  gaa  conatant,  w  ia  the  angu¬ 
lar  frequency  of  the  acouatic  wave,  and  t  repreaenta  time. 

Entropy  verlatione  are  detected  aa  temperature  fluctuatlona  given  by 
the  relation 


T(x,t)  -  T _{l  +  2e  ain  -j-  coa  (mt  -  (5.4) 

8  8 

where  Tfia  the  flame  temperature,  x  ia  the  distance  from  the  propellant 
aurface,  and  Vg  ia  the  gas  velocity  (Ref.  37) 

It  waa  noted  in  Ref.  37  that  at  any  lnatant  in  time  there  would  not 
only  be  a  ainuaoldal  temperature  variation  in  apace,  with  amplitude 
Tf  c  (y-1)/y,  but  that  there  would  alao  be  a  time-varying  temperature 
observed  at  fixed  pointa  on  the  burner.  In  addition,  if  one  were  to  ob- 
eerve  the  magnitude  of  temperature  fluctuatlona  at  varloua  pointa  along 
the  burner,  it  would  be  found  that  the  amplitude  of  temperature  fluctu- 
ationa  variea  with  dlatance  from  the  propellant  aurface.  Maxima  in  the 
temperature  fluctuation  would  be  found  at  locations  given  by 


V  ■  (n  +  w  vgr’ 

max  6 


(5.5) 


while  the  temperature  fluctuaticn  would  be  zero  at  pointa  located  at 


x_ i  -  n  v  T,  (5.6) 

min  8 

where  n  repreaenta  any  poaltive  integer  and  T  ia  the  period  of  acouatic 
oscillation.  The  amplitude  of  temperature  fluctuation  (Eq.  5.5)  la 
equal  to  2  Tf  c(y-l)/y,  i.e.,  twice  the  value  of  the  spatial  temperature 
variation. 

Equations  5.4,  5.5,  and  5.6  provide  a  standard  for  comparing  the 
foregoing  theory  with  experimental  obaervationa.  The  observation  of  gaa 
temperature  variations  aa  a  function  of  apace  and/or  time  coupled  with 
knowledge  of  the  acouatic  frequency  and  the  gaa  velocity  ahould  be  suffi¬ 
cient  to  verify  exlatence  of  the  entropy  wave  phenomenon.  Sufficiently 
detailed  knowledge  of  the  relationship  between  the  temperature  fluctua¬ 
tions  and  the  acoustic  pressure  would  also  permit  a  determination  of  the 
o/c  term  through  the  "temperature  response  function"  aa  diacuased  in 
Ref.  37  and  38. 

Several  attempts  have  been  made  in  the  past  to  detect  nonisentropic 
behavior  but  theae  efforts  did  not  yield  a  positive  result.  Recent 
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theoretical  studies  Indicated  that  the  previous  experiments  may  have 
failed  because  they  were  conducted  at  frequencies  111-sulted  for  display¬ 
ing  of  the  phenomenon  (Ref.  39).  Since  the  NWC  5.5-in.  burner  was  capa¬ 
ble  of  operating  within  the  frequency  range  of  Interest,  it  was  decided 
to  attempt  to  obtain  evidence  of  nonlsentroplc  behavior  In  that  apparatus. 

Two  types  of  experiments  were  conducted:  (1)  a  series  In  which 
photographic  recording  was  utilized  to  obtain  space-tine  records  of  the 
light  intensity  variations  In  ths  gas  after  It  left  the  combustion  zone; 

(2)  a  series  in  which  high-speed  temperature  measurements  was  made 
using  an  adaptation  of  the  brlghtness-emisslvity  technique.  Three  pro¬ 
pellant  compositions  were  used  in  the  two  seriee  of  experiments:  (1)  a 
bimodal  polyurethane-AP  composition  prepared  at  Princeton  University 
(deslgneted  PU  in  this  report),  and  cast  into  metal  cups,  (2)  a  modified 
version  of  A-35  containing  0.4X  KaCl  which  was  designated  A-152,  and 

(3)  A-13  propellant.  Comoosltlon  of  these  propellants  appears  In  Table 
5.1. 

5.3.2.  Photographic  Studies 

The  series  of  photographic  experiments  utilized  a  special  burner 
test  section  which  was  fitted  with  a  nitrogen-flushed  slit  window  (Fig. 
5.5).  An  optical  bench  was  used  to  mount  a  70  mm  streak  camera  and  a 
long  focal  length  objective  lens.  A  cathode  ray  oscilloscope  and  mirror 
were  arranged  so  that  the  oscilloscope  ecreen  was  focused  on  the  film 
plane  in  the  camera  In  line  with  the  image  of  the  window  slit  image. 

Thus  the  acoustic  pressure  trace,  which  appeared  on  the  cathode  ray 
screen,  was  recorded  simultaneously  with  the  radiant  intensity  variations 
that  appeared  in  the  slit  window  on  the  burner.  This  arrangement  per¬ 
mitted  a  detailed  comparison  to  be  made  between  events  occurring  in  the 
burner  gas  and  the  acoustic  pressure.  A  block  diagram  of  the  equipment 
arrangement  used  in  the  photographic  studies  is  shown  In  Fig.  5.6. 

A  portion  of  streak  camera  record  is  shown  in  Fig.  5.7.  The  acous¬ 
tic  pressure  and  approximate  poeition  of  the  propellant  surface  have 
been  added  to  the  figure  as  an  aid  in  interpretation.  Major  features  of 
interest  have  been  inked  onto  the  figure  to  the  left  of  the  photographic 
record  for  the  same  reason.  The  photographic  portion  of  this  figure  is 
a  positive  ao  that  high  and  low  temperature  regions  in  the  figure  appear 
as  light  and  dark  areas  respectively.  Calculations  from  Eq.  5.5  and  5.6, 
based  on  an  estimate  of  the  mean  value  of  the  gas  velocity  as  it  leaves 
the  propellant  surface,  were  made  in  order  to  locate  the  positions  of 
the  maximum  and  minimum  temperature  variatione.  These  are  deslgneted  in 
the  figure  as  T’m  x  and  T’m^n  respectively.  Note  also  that  a  time  and 
distance  scale  are  provided  in  Fig.  5.7  and  that  time  runs  from  right  to 
left. 

Major  features  in  the  reconstructed  portion  of  the  record  are  iden¬ 
tified  with  letters  as  follows: 
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FIG.  5.5.  Test  Section  Containing  Nitrogen  Flushed 
Slit  Window  Used  on  the  5.5-Inch  T-Burner  for  Photo¬ 
graphic  Studies  of  Nonlsentroplc  Phenomena  In  Low 
Frequency  Acoustic  Instability. 


FIG.  5.6.  Blo;k  Diagram  of  Streak  Carnets  Arrange¬ 
ment  Used  for  photographic  Study  of  Nonlsentroplc 
Phenomena  In  the  5.5-Inch  Diameter  Burner. 


A.  Periodic  vertical  darkened  zones  appear  In  the  record  which 
occur  st  times  when  the  acoustic  pressure  Is  a  minimum.  These  areas 
appear  to  be  related  to  adiabatic  cooling  of  gas  throughout  the  burner 
section  covered  by  the  window  caused  by  lowering  of  the  acoustic  pres¬ 
sure.  (Note  that  there  Is  also  a  general  brightening  at  acoustic  pres¬ 
sure  maxima  related  to  adiabatic  hestlng  from  compression  of  the  gas  In 
the  burner  by  the  acoustic  wave.) 

B.  Dark  diagonal  lines,  designated  "B",  are  relatively  persistent 
In  their  travel  acrocs  the  section  covered  by  the  window.  These  repre¬ 
sent  elements  of  relatively  cold  gas  which  appear  to  have  originated 
from  the  propellant  surface  st  a  time  when  the  acoustic  pressure  was  a 
minimum.  This  assertion  of  the  origin  of  theie  elements  of  gas  Is 
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FIG.  5.7.  Portion  of  Streak  Camera  Record  from  Teat  with  PU  Pro¬ 
pellant.  Mean  burner  preaaure  was  75  paia,  acoustic  pressure 
amplitude  waa  9  psi  p-p,  time  after  Ignition  was  2.3  see.  The 
position  of  the  burning  surface,  the  acoustic  preaaure,  and  some 
of  the  major  features  of  the  record  have  been  reconstructed  for 
clarity.  Points  of  particular  interest  in  the  reconstructed 
portion  of  the  record  ere  designated  by  lettera  in  the  figure, 
with  appropriate  description  in  the  test.  Locations  of  maximum 
and  minimum  temperature  fluctuations,  designated  T’,^  and  T'nin 
respectively,  were  calculated  from  Eq.  5.5  and  5.6. 


based  on  an  extrapolation  of  the  information  on  the  photographic  record 
back  to  the  propellant  surfsce  as  indicated  in  the  inked  reconstruction 
to  the  left  of  the  photographic  portion  of  the  figure.  It  waa  assumed 
in  making  the  extrapolation  that  the  gas  velocity  in  the  region  of  ex¬ 
trapolation  is  relatively  unaffected  by  the  acoustic  wave.  If  the  in¬ 
terpretation  given  here  ia  correct,  a  phenomenon  is  present  in  these 
tests  that  the  theory  in  Section  5.3.1  does  not  take  into  account. 

C.  Apparent  inflection  points  occur  in  the  motion  of  cold  gas 
elements  (designated  "C".  The  photographic  record  auggests  that  a  aig- 
nif leant  fluctuation  occurs  involving,  perhaps,  a  momentary  flow  re¬ 
versal.  It  should  be  noted  that  the  theory  described  in  Section  5.3.1 
assumes  uniform  gaa  flow  and  does  not  provide  for  a  flow  perturbation  of 
the  sort  indicated  in  the  record. 

D.  Dark  areas  designated  "D"  represent  regions  of  cold  gas,  which 
at  any  inatent  in  time  occupy  a  position  midway  between  "8"  elements  of 
gas.  Since  these  appear  between  areas  of  gas  ("B"  designation)  which 
originated  at  acoustic  pressure  minima,  it  seems  that  the  "D"  regions 
sre  related  to  the  entropy  fluctuationa  referred  to  in  Section  5.3.1. 


76 


NWC  TP  4478 


There  are  several  difficulties  related  to  the  above  discussion. 

First  is  the  problem  of  identification  of  the  cold  gas  that  appears  to 
leave  the  propellant  surface  at  an  acoustic  pressure  minimum.  The  ex¬ 
trapolation  of  an  element  of  gas  back  to  the  surface  has  to  be  done  over 
more  than  an  inch  for  which  there  is  no  photographic  record.  In  doing 
the  extrapolation,  the  assumption  has  been  made  that  the  acoustic  velo¬ 
city  perturbation  in  this  region  is  negligible  since  it  is  so  close  to 
the  propellant  surface,  relative  to  the  acoustic  wavelength.  Second, 
the  film  was  overexposed  causing  a  diminishing  of  the  relative  import¬ 
ance  of  cooler  regions  in  the  gss.  Some  of  the  effects  of  overexposure 
have  been  counteracted  by  printing  on  high-contrast  material  but  this 
may  have  introduced  further  distortion  of  the  quality  of  the  photograph. 
Problems  in  photographic  record  interpretation  can  probably  be  resolved 
by  repeating  tests  with  altered  burner  parts  which  allow  the  propellant 
to  be  moved  up  to  the  end  of  the  window  and  by  using  a  camera  lens  with 
smaller  aperture  than  was  used  in  the  first  test  series. 

5.3.3.  High-Speed  Temperature  Measurements 

The  second  part  of  this  investigation  involved  acquisition  of  high¬ 
speed  gas  temperature  data  from  a  region  near  the  propellant  surface. 

The  measurement  technique  utilized  an  adaptation  of  the  brightness- 
emissivity  approach  mentioned  in  Ref.  37.  The  apparatus  is  shown  in 
Fig.  5.8  and  consists  of:  (A)  a  ribbon  filament  lamp  of  known  bright¬ 
ness  temperature,  (B)  a  syatem  of  lenses  and  slits  to  focus  and  colli¬ 
mate  the  light,  (C)  a  rotating  perforated  metal  wheel  which  periodically 
interrupts  light  from  the  known  source,  (D)  a  special  test  section  which 
fits  between  the  burner  and  an  end  cap  containing  propellant  that  con¬ 
tains  two  small  nitrogen-flushed  windows  which  permit  light  from  the 
known  source  to  pass  through  the  burner,  and  (E)  a  monochromater 
equipped  with  a  photomultiplier  tube  which  receives  and  detects  radi¬ 
ation  in  the  sodium  D-line  portion  of  the  optical  spectrum.  The  photo¬ 
multiplier  signal  was  recorded  on  an  oscillograph  simultaneously  with 
the  chamber  pressure.  Prior  to  each  run,  the  source  was  calibrated  with 
a  precision  pyrometer  and  a  short  calibration  record  was  obtained  on  the 
oscillograph  giving  the  deflection  due  to  the  lamp  intensity.  During  the 
test,  the  radiation  received  by  the  photomultiplier  tube  alternated  be¬ 
tween  that  from  the  flame  and  that  of  the  flame  plus  radiation  from  the 
source.  A  portion  of  a  test  record  is  shown  in  Fig.  5.9. 

Data  reduction  following  a  test  was  accomplished  by  making  deflec¬ 
tion  measurements  of  the  pre-run  source  calibration,  and  measuring  the 
signal  deflections  from  the  photomultiplier  and  the  pressure  transducer 
over  one  or  more  cycles  of  oscillation  at  selected  portions  of  the  test. 

A  typical  plot  of  temperature-pressure  data  is  shown  in  Fig.  5.10. 

Data  from  selected  tests  appear  in  Table  5.2.  Columns  {a]  through 
[ h]  represent  data  derived  from  the  oscillograph  recordings.  Alphabeti¬ 
cal  suffixes  to  the  test  number  are  used  to  identify  the  portion  of  the 
test  record  from  which  the  data  was  taken.  The  time  at  which  the  data 
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FIG.  5.8.  Apparatus  Used  for  Brlghtness-Emlsslvlty  Temperature 
Determination.  Lettered  Items  are  described  In  the  text. 
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FIG.  5.9.  Portion  of  Test  Record  From  a  Test  in 
Which  High-Speed  Temperature  Data  was  Acquired. 


FIG.  5.10.  Temperature-Pressure  Data  Reduced 
and  Plotted  on  a  Common  Time  Base  From  Test 
Record  Similar  to  That  Shown  in  Fig.  5.9. 
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was  taken  appears  in  column  [a]  and  is  the  time  elapsed  after  ignition. 
Frequency  was  determined  by  measurement  of  the  oscillograph  record 
(column  [b])  while  the  pressure  and  temperature  data  were  obtained  from 
plots  similar  to  that  shown  in  Fig.  5.10.  Column  [gj  was  calculated  by 
taking  the  ratio  of  the  pressure  oscillation  amplitude  to  the  mean 
pressure.  The  phase  between  temperature  and  pressure  (column  [h])  waa 
computed  from  measurements  of  the  pressure-temperature  plots  (as  in 
Fig.  5.10).  Columns  [i]  and  [j]  are  computed  temperature  fluctuations 
thst  will  be  discussed  in  some  detail  below. 

In  general,  the  tabulsted  pressure-temperature  data  show  quite  a 
diversity  of  behavior.  The  A-152  data  represent  the  most  comprehensive 
complilstion  of  information  available  at  present,  but  the  question  of 
whether  the  data  fits  a  consistent  pattern  must  await  completion  of  the 
data  analysis.  One  outstanding  feature,  for  which  there  is  no  explana¬ 
tion  at  present,  is  the  surprisingly  high  amplitude  of  temperature  os¬ 
cillation  observed  in  the  testa  above  40  cps.  A  comparison  of  observed 
temperature  fluctuation  amplitude  with  calculated  amplitudes,  baaed  on 
two  rather  different  assumptions  suggests  that  neither  of  the  explana¬ 
tions  for  the  temperature  oscillations  advanced  so  far  can  explain  what 
is  observed  in  the  burner.  The  first  comparison  is  with  the  temperature 
fluctuation  predicted  on  the  assumption  thst  the  gas  in  the  burner  be¬ 
haves  isentropically.  This  assumption  leads  to  the  calculated  temper¬ 
atures  entered  in  column  [i].  Comparison  of  the  lsentroplc  calculation 
and  the  observed  temperature  oscillations  indicates  the  increased 
difference  between  the  lsentroplc  prediction  and  the  observed  vslues. 

If  the  assumptions  of  Section  5.3.1  were  valid,  the  observed  tempera¬ 
ture  oscillation  amplitude  would  not  exceed  the  lsentroplc  value  by 
more  than  a  fsctor  of  two.  As  one  can  aee,  the  observed  amplitudes 
from  tests  with  frequencies  above  40  cpa  far  exceed  the  expected  range. 

Another  fscet  of  the  temperature  oscillations  was  considered, 
namely  the  effect  of  a  fluctuating  flame  temperature.  It  was  noted, 
in  comparing  mean  temperature  and  mean  pressure  from  PU  snd  A-152  data, 
that  the  temperature  was  a  function  of  pressure.  The  dependence  of 
temperature  on  pressure  for  these  two  propellants  is  indicated  in 
Fig.  5.11.  Assuming  that  the  temperature  follows  the  pressure  oscil¬ 
lations  in  the  ssme  manner  as  indicated  in  Fig.  5.11,  one  can  aaslgn 
a  sensitivity  of  temperature  to  pressure  which  can  be  used  to  compute 
the  tempersture  oscillation  caused  by  an  observed  pressure  fluctuation. 
Computations  based  on  the  above  assumptions  are  entered  in  column  [j] 
of  Table  5.2  and  the  temperature  sensitivity  for  the  two  propellants 
appears  in  a  footnote  in  that  table.  Comparison  of  the  observed  and 
computed  dsts  indicates  a  large  disparity  and  that  the  temperature 
sensitivity  to  pressure  does  not  offer  an  adequate  explanation  of  the 
large  temperature  oscillations  observed  at  the  higher  frequencies. 

A  detailed  report  is  in  preparation  which  describes  the  photo¬ 
graphic  and  hlgh-apeed  temperature  measurements  in  greater  detail. 

The  in  formation  will  be  released  as  a  NWC  report  when  completed. 
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FIG.  5.11.  Plot  of  Measured  Mean  Gas  Temperature 
end  Mesn  Pressure  for  PU  and  A-152  Propellants. 


5.3.4.  A  New  Mathematical  Model 

Detailed  Interpretation  of  the  streak  photography  and  the  temper¬ 
ature  data  Is  hampered  by  presence  of  phenomena  not  sccounted  for  In 
existing  theories.  It  was  felt  that  a  mathematical  model  could  be  de¬ 
vised  which  would  account  for  a  variable  flame  temperature  end  non- 
uniform  flow  of  gsses.  Steps  In  that  direction  were  taken  lste  in  the 
report  period. 

A  new  theory  is  presented  here  which  describes  the  history  of  a 
gas  particle  from  the  time  it  leaves  the  reaction  zone  as  a  function  of 
a  parameter  ♦  which  denotes  the  phasing  of  the  pressure  and  velocity 
fields  at  the  emergence  of  the  particle  from  the  combustion  zone. 

The  oscillating  components  of  pressure  and  velocity  were  treated 
as  first  order  perturbation  quantities  and  since  the  amplitude  of  the 
oscillating  component  of  velocity  wss  of  the  order  of  the  mean  velocity, 
a  point  clearly  shown  by  the  evidence  of  flow  reversal  In  some  of  the 
photographic  test  records,  the  nonconvectlve  wave  equation  was  used  to 
describe  the  acoustic  fields.  Pressure  antinodes  at  each  end  of  the 
chamber  were  used  as  boundary  conditions  with  the  resulting  oscillating 
pressure  and  velocity  fields  given  by 


p’  ■  p  c  cos  (sx/L)  cos  (wt  +  ♦)  (5.5) 
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and 


v'  ■  V  sin  (itx/L)  sin  (wt  +  i|>)  (5.6) 

respectively.  In  Eq.  5.5  and  5.6  p  is  the  mean  pressure,  c  is  the  ratio 
of  the  maximum  value  of  the  pressure  perturbation  to  mean  pressure,  V  is 
the  velocity  amplitude,  x  is  the  distance  from  the  reaction  zone  snd  L  is 
the  T-burner  length.  The  particle  velocity  is  given  by 


v  «  v  +  v' 


(5.7) 


and  the  relationship  between  temperature  and  pressure  given  by 


T 

T 


r 


(5.8) 


in  which  pr  was  the  pressure  at  a  reference  plane  immediately  outside  of 
the  flame  when  the  particle  in  question  entered  the  acoustic  field  and 
Tr  is  the  temperature  at  the  reference  plane,  which  is  assumed  to  be 
constant  for  this  discussion.  The  reference  pressure  given  by 


p  “  P  +  pc  cos  ^ 


(5.9) 


was  obtained  by  adding  the  mean  pressure  p  to  the  acoustic  pressure  with 
x  and  t  set  equal  to  zero.  The  following  expression  for  relative  tem¬ 
perature 


T’  +  T 

_ r 

T 

r 


.cos  (flx/L)  cos  (ut  +  i|>)  -  cos  <i. 
1  +  e  cos 


.(5.10) 


was  obtained  by  substitution  of  Eq.  5.7  snd  5.9  into  Eq.  5.8.  In  order 
to  evaluate  T/Tr  it  is  necessary  to  determine  the  relationship  between 
x  and  t  for  s  given  if>.  This  is  readily  obtained  from  the  relation 


x 


I"1 


t 

|  [v  +  V  sin  sin  (at  +  ♦)]  dt 

o 


(5.11) 
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Solution  of  Eq.  5.11  thus  gives  the  trajectory  of  a  gas  particle  leaving 
the  aurface  at  tine  zero  for  the  assumed  boundary  conditions  and  phase 
angle  i|i.  Equations  5.10  and  5.11  together  give  the  temperature  hlatory 
and  the  trajectory  of  a  gas  particle  subsequent  to  leaving  the  surface 
as  a  function  of  the  parameter 

Additional  steps  In  modifying  the  analysis  to  provide  a  more  ade¬ 
quate  description  of  significant  parameters  In  the  problem  of  nonlsen- 
troplc  behavior  need  to  be  undertaken.  These  Include  replacing  the 
rigid  solid  boundary  condition  at  the  combustion  zone  with  an  scouetlc 
admittance  and  replacing  the  constant  flame  temperature  with  one  which 
la  pressure  sensitive. 

It  la  planned  to  present  solutions  of  the  equations  in  the  form  of 
dimensionless  plots.  With  these  available,  particle  trajectories  and 
constant  temperature  lines  may  then  be  superimposed  on  experimental  test 
records,  t^e  phase  angle  being  readily  obtained  from  the  pressure 
record.  The  relative  brightness  patterns  In  the  streak  data  will  then 
provide  a  qualitative  test  of  the  mathematical  model. 


5.4.  SUMMARY 

Tests  to  determine  range  of  unstable  combustion  behavior  on  one 
metallized  and  two  unmetsllized  propellants  wse  accomplished  prior  to 
more  definitive  testing.  It  is  anticipated  that  extensive  responee 
function  testing  on  the  two  unmetallized  compositions  will  be  accomplished 
in  the  near  future. 

Preliminary  testing  with  a  preseed  AP-A1  composition  in  the  5.5- 
lnch  diameter  T-burner  demonstrated  the  feasibility  of  using  such  mate- 
riale  in  lnveetigatlng  the  role  of  metal  combustion  In  preferred- 
frequency  Instability. 

Tnveetlgatlon  of  nonleentroplc  behavior  In  the  low  frequency  range 
revealed  that,  under  the  conditions  teeted,  the  flame  temperature  may 
fluctuate  elgniflcsntly  with  the  acouetlc  preaeure.  While  It  appeared 
that  the  flame  temperature  of  two  of  the  propellants  tested  wae  a  func¬ 
tion  of  the  mean  preeeure,  the  resulting  temperature  variation  wae  not 
large  enough  In  many  caeee  to  explain  the  observed  temperature  fluctu¬ 
ations.  The  magnitude  of  the  temperature  oscillations  and  the  evidence 
from  the  photographic  recorde  atrongly  auggeeted  that  factore  other 
than  thoee  dlscueeed  above  were  Involved.  One  plaueible,  but  ae  yet 
unevaluated,  explanation  wae  that  the  binder  and  oxidizer  were  consumed 
at  different  ratee  which  vary  over  a  cycle  of  preesure  oecillatlon, 
reeulting  in  a  fluctuation  in  the  fuel-to-oxidizer  ratio. 

The  observatlone  of  nonleentroplc  behavior  etrongly  suggeeted  that 
preeent  mathematical  models  which  eeek  to  describe  unstable  combustion 
may  need  additional  reflnemente  in  order  to  accuretely  deecribe  the 
phenomena  that  appzar  to  be  Involved. 
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6.  SUMMARY 


Each  of  the  preceding  sections  contains  its  own  extensive  discus¬ 
sion,  and  the  major  accomplishments  of  each  section  are  summarized  here 
in  order  to  provide  a  comprehensive  view  of  the  entire  program. 


6.1.  DECOMPOSITION  AND  DEFLAGRATION  OF  AMMONIUM  PERCHLORATE 

Past  and  current  work  has  shown  that  the  physical  character  of  the 
surface  residue  of  AP  is  significantly  dependent  on  the  imposed  heating 
rate.  Ammonium  perchlorate  apparently  sublimes  under  conditions  of  iso¬ 
thermal  decomposition  and  there  is  no  evidence  of  a  surface  melt;  under 
conditions  of  rapid  heating  such  as  burning  in  a  methane  atmosphere,  a 
surface  melt  is  observed  on  the  AP;  and  in  the  case  of  self-sustaining 
deflagration,  the  surface  melt  is  frothy  due  to  the  entrapment  of  tran¬ 
spiring  gases.  This  frothy  melt  is  probably  the  site  of  the  major  heat 
release  necessary  to  sustain  deflagration. 


6.2.  BENAVIOR  OF  ALUMINUM  IN  COMPOSITE  PROPELLANTS 

Current  work  has  indicated  that  an  Important  step  in  the  aluminum 
agglomeration  (sintering  together)  process  is  characterized  by  the  frac¬ 
ture  of  the  aluminum  oxide  coating  (below  its  melting  point)  due  to 
stress  produced  by  the  differential  thermal  expansion  of  aluminum  and 
aluminum  oxide  as  the  particles  are  heated  by  the  combustion  wave  and 
oxidative  self-heating.  Molten  aluminum  flows  through  the  cracka,  joins 
with  adjacent  particles  and  is  oxidized,  thus  forming  aluminum/aluminum 
oxide  particles.  By  including  this  Important  step  of  the  agglomeration 
process,  it  is  possible  to  give  a  detailed,  qualitative  description  of 
the  behavior  of  aluminum  in  composite  propellants  from  the  unheated  con¬ 
dition  in  the  propellant  to  the  condition  of  molten  balls  of  aluminum 
with  oxide  "caps"  burning  by  a  vapor  transport  mechanism  in  the  gas 
phase  above  the  propellant  surface. 


6.3.  NONACOUSTIC  COMBUSTION  INSTABILITY 

Two  analytical  modela  of  solid  propellant  combustion  instability 
applicable  to  the  low  frequency  regime  have  been  reduced  to  a  common 
notation,  were  found  to  be  numerically  identical,  and  have  been  compared 
with  experimental  data.  Only  qualitative  agreement  between  experimental 
data  and  model  predictions  was  found,  and  even  then  several  large 
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discrepancies  were  discovered.  One  discrepancy  associated  with  the 
heterogeneous  nsture  of  the  propellant  was  experimentally  explored  in 
detail.  A  simple,  phenomenological  argument  relating  oscillation  fre¬ 
quency  to  oxidizer  particle  size  was  developed  and,  when  applied  in  con¬ 
junction  with  the  analytical  models,  was  found  to  be  useful  in  inter¬ 
preting  the  experimental  data. 


6.4.  ACOUSTIC  INSTABILITY 

Experiments  with  pressed  AP-A1  material  in  the  NWC  5.5-inch  T-burner 
indicate  that  auch  material  can  be  successfully  used  to  investigate  the 
interaction  of  acoustic  waves  with  spontaneously  occurring  combustion  os- 
cillstions. 

Investigation  of  nonlsentropic  behavior  has  shown  that  pehnomena 
were  Involved  which  were  not  considered  in  presently  svsllsble  mathe¬ 
matical  models.  Experimental  information  acquired  during  this  report 
period  indicated  that  flame  temperature  variations  and  gas  velocity  per¬ 
turbations  were  of  the  same  order  of  importance  as  the  temperature  fluc¬ 
tuations  Induced  by  variations  in  entropy.  Steps  were  initiated  to  syn¬ 
thesize  a  more  relevant  mathematical  model  which  will  include  consider¬ 
ation  of  nonuniform  gaa  velocity  and  a  variable  flame  temperature,  as 
well  as  Improvements  in  the  methods  of  obtaining  experimental  measure¬ 
ments. 
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Appendix  A 

BURNI..J  SATE  DATA  AND  TABULAR  EXPERIMENTAL  DATA  OBTAINED 
WITH  RESPECT  TO  NONACOUSTIC  COMBUSTION  INSTABILITY 


TABLF  A-l.  Data  Obtained  WilS  A-146  Propellont. 


Run 

nq. 

L*.  in 

Frequency, 

cp> 

Prvtiurv, 

ptio 

Burning 

rot* 

in/sec 

Growth 

con»tor»t, 

TCf.4 

meec 

DimonsionUs* 

frtqutncy 

afa> 

41TC 

Phq*e 

ihift 

rodiqni 

Root  pa 
of  the 

reepons 

functiq 

906 

20 

168 

168 

.162 

.850 

7.3 

.897 

32 

55 

195 

.168 

1.360 

2.2 

.469 

907 

18 

162 

145 

.154 

. 

.765 

7.8 

.777 

• 

64.8 

158 

.159 

1.904 

. 

908 

15.9 

165 

141 

.153 

97.3 

.676 

8.0 

.699 

.580 

1.07 

17.1 

150 

148 

.155 

114 

.727 

7.1 

.684 

.563 

1.08 

19 

131 

153 

.157 

39.1 

.808 

6.1 

.661 

.570 

1.03 

21 

131 

153 

.157 

-37.8 

.893 

6.1 

.731 

.647 

.936 

910 

14.6 

172 

121 

.144 

113 

.621 

9.5 

.669 

.559 

1.07 

16.1 

150 

130 

.148 

54.3 

.684 

8.3 

.679 

.578 

1.04 

18.5 

133 

131 

.149 

36.7 

.786 

6.8 

.655 

.567 

1.03 

24.2 

104 

130 

.148 

. 

1.029 

5.4 

.671 

899 

45 

56 

95 

.129 

27.2 

1.913 

3.8 

.671 

.568 

1.05 

901 

13.8 

171 

115 

.141 

.587 

?.8 

.629 

15.1 

165 

125 

.1465 

55.5 

.642 

8.8 

.665 

.571 

1.04 

16.8 

133 

125 

.1465 

70 

.714 

7.1 

.594 

.515 

1.05 

911 

14.0 

158 

120 

.143 

80.5 

.595 

8.9 

.590 

.513 

1.05 

15.1 

148 

125 

.1465 

85 

.642 

7.9 

.595 

.513 

1.05 

16.5 

139 

127 

.147 

61.8 

.701 

7.3 

.612 

.531 

1.04 

18.9 

127 

124 

.146 

38.5 

.803 

6.8 

.642 

.557 

1.03 

23.5 

105 

121 

.144 

...... 

.999 

5.8 

.659 

28.6 

97 

122 

.145 

1.216 

5.3 

.739 

916 

12.4 

160 

98 

.131 

139 

.527 

10.7 

.529 

.458 

1.07 

13.0 

162 

105 

.135 

126 

.553 

10.2 

.562 

.483 

1.07 

915 

14.1 

148 

102 

.133 

95 

.599 

9.6 

.555 

.483 

1.06 

15.8 

128 

106 

.136 

96 

.672 

7.9 

.539 

.468 

1.07 

916 

23.6 

100 

93 

.127 

67 

1.003 

7.1 

.628 

.532 

1.07 

26.2 

98 

101 

.133 

53.3 

1.114 

6.3 

.683 

.573 

1.06 

917 

29.3 

83 

98 

.131 

*45.8 

1.245 

5.5 

.648 

.550 

1.06 

918 

39.2 

72 

96 

.130 

-17.5 

1.666 

4.875 

.753 

.660 

0.97 

42.4 

58 

85 

.122 

-  5.3 

1.802 

4.450 

.655 

.584 

1.00 
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TA8LE  A-2.  Data  Obtained  With  A-148  Propellant. 


Run 

no. 

L*,  in 

Frequency, 

cp* 

ptio 

Burning 

rote 

in/tec 

Growth 
con  ktont, 

rch. 

mkee 

Dimen  >  ion  lets 
frequency 

a.fb) 

T~ 

d)Tc 

Reol  port 
of  the 
rokponk# 
function 

1001 

19 

150 

114 

.116 

.806 

13.2 

.759 

. 

20.1 

125 

114 

.116 

.855 

11.0 

.671 

127 

115 

.117 

66.2 

.98 

11.2 

.781 

1.07 

23.8 

119 

113 

.116 

57.2 

1.01 

10.5 

.754 

1.06 

24.9 

116 

111 

.115 

36.0 

1.06 

10.4 

.772 

1.04 

27.2 

108 

108 

.115 

22.3 

1.16 

9.76 

.786 

1.03 

31.3 

100 

105 

.113 

0 

1.33 

9.23 

.835 

1.00 

32.7 

98 

104 

.113 

0 

1.39 

9.04 

.855 

1.00 

40.2 

73 

90 

.107 

28.7 

1.71 

7.56 

.783 

1,05 

1002 

25 

100 

81 

.102 

. 

1.06 

11.4 

.666 

25.6 

~I00 

85 

.105 

1.09 

11.0 

.686 

26.5 

97 

85 

.105 

...... 

1.13 

10.7 

.688 

27.1 

67 

82 

.103 

1.15 

9.39 

.628 

26.2 

91 

84.5 

.105 

1.20 

9.99 

.685 

...... 

26.5 

89 

84 

.104 

1.26 

9.78 

.704 

30.6 

85 

86 

.106 

1.30 

8.70 

.694 

32.2 

60 

83 

.103 

1.37 

8.48 

.688 

. 

33.6 

81 

80 

.102 

1.43 

9.47 

.728 

35.0 

60 

60 

.102 

1.49 

9.34 

.748 

1003 

33 

72 

70 

.0955 

1.42 

9.49 

.642 

. 

36.6 

70 

.0955 

1.56 

40.2 

65 

.0915 

34.3 

1.71 

1.06 

45 

57 

67 

.0935 

27.2 

1.91 

7.70 

.684 

1.05 

1124 

17.5 

115 

91 

.107 

■ 

.74 

11.5 

.534 

16.2 

120 

94 

.108 

78 

11.8 

.588 

19 

110 

93 

.108 

.81 

10.8 

.560 

19.7 

120 

93 

.108 

.83 

11.8 

.626 

20.7 

94 

.108 

.88 

22 

94 

.106 

.94 

1125 

28.1 

80 

62 

.104 

■ 

1.20 

8.48 

.602 

29.8 

80 

94 

.108 

84.0 

1.27 

7.83 

.638 

1.11 

33.8 

87 

97 

.109 

29.2 

1.43 

8.35 

.781 

1.04 

40.6 

75 

94 

.106 

25.6 

1.73 

7.35 

.815 

1.05 

1126 

44.7 

71 

94 

.106 

19.0 

1.90 

6.96 

.847 

1.04 

1127 

19 

100 

103.5 

.113 

28.2 

.809 

8.92 

.508 

1.02 

22 

no 

107.5 

.113 

5.2 

.94 

9.81 

.650 

1.01 

25.9 

105 

105.0 

.112 

-16.2 

1.10 

9.56 

.725 

.98 
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18.3 

115 

107 

.115 

.78 

9.69 

.563 

20.4 

133 

117 

.117 

76.7 

.87 

11.1 

.726 

1.07 

25.3 

~  6.7 

120 

.118 

0 

i.o; 

.55 

.045 

1.00 

1129 

50.7 

72 

93.5 

.110 

2.16 

6.78 

.976 
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TABLE  A-3.  Onto  Obtained  With  A- 149  Propellant. 


Run 

no. 

L\  In 

Fr*qu*ncy, 

cps 

Rr*ssur*, 

psio 

Burn  in  p 
rot* 

in/s*c 

Growth 

constant, 

s*c-l 

rc  h 
ms*c 

Dim*nsiont*ss 

fr*qu*ncy 

df<U 

airc 

Roof  port 
of  th* 
r* spans* 
function 

TJ 

1016 

106.3 

10.8 

68 

.064 

.... 

4.51 

3.0 

.306 

116 

10.7 

69 

.064 

4.93 

3.0 

.332 

1019 

41 

130 

.084 

...... 

1.74 

4.86 

.326 

...... 

54.5 

150 

.089 

2.32 

3.38 

.341 

...... 

61.4 

158 

.091 

2.61 

2.77 

.328 

64.4 

27.3 

158 

.091 

..... 

2.74 

3.78 

.470 

1020 

57.5- 

. 

151.5 

.089 

2.42- 

..... 

107 

4.52 

1021 

35 

18.8 

108 

.077 

..... 

1.49 

.176 

...... 

37.7 

25 

110 

.078 

1.60 

HI 

.251 

— 

1022 

60 

20 

103 

.075 

2.55 

.321 

1023 

70.3 

13.3 

55 

.058 

5.42 

2.99 

1 

.249 

1.04 

69 

13.2 

58 

.059 

-4.46 

3.78 

.314 

.98 
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TABLE  A-7.  Doto  Obtained  With  A-157  Propellant. 


Run 

no. 

L*,  in 

Frequency, 

cps 

Pressure, 

psio 

Burning 

rote 

in/eec 

Growth 

constant, 

sec-' 

Tzh. 

msec 

Dimensionless 

frequency 

a-tCLi 

aire 

Real  part 
of  the 
response 
function 

73 

1159 

31.9 

46 

76 

.117 

27.3 

1.36 

3.84 

.393 

1.037 

1151 

53 

33 

69.5 

.110 

2.25 

3.12 

.466 

-... 

1163 

24.1 

62.5 

73.5 

.115 

1.02 

5.40 

.401 

..... 

25.9 

70.0 

83.5 

.120 

1.10 

5.56 

.484 

1194 

30 

47 

93.5 

.120 

31.5 

1.28 

3.73 

.378 

1.04 

1197 

27.1 

51.5 

88.5 

.120 

31.8 

1.15 

4.09 

.372 

1.04 

30.0 

55.0 

88.5 

.120 

32.0 

1.28 

4.36 

.443 

1.04 

31.6 

61 

90.5 

.120 

1.34 

4.85 

.514 

..... 

34.9 

50 

87.0 

.120 

38.4 

1.48 

3.97 

.465 

1.06 

39.0 

so 

83.5 

.120 

37.0 

1.66 

3.97 

.523 

1.06 

1198 

45.6 

22.8 

63.0 

.103 

1.94 

2.46 

.278 

— 

1700 

68 

27.5 

63.5 

.103 

2.89 

2.97 

.499 

73 

33.3 

60.5 

.101 

3.10 

3.72 

.649 

— 

1202 

66.1 

36.0 

72.5 

.114 

23.4 

2.81 

3.16 

.635 

1.07 

74.1 

35.2 

75.5 

.115 

3.15 

3.04 

.696 

79.3 

31.0 

58.5 

.095 

— 

3.37 

3.94 

.656 

1211 

20.1 

85 

98.5 

.120 

...... 

.85 

6.75 

.454 

. 

27.0 

75 

96 

.120 

1.15 

5.96 

.542 

. 

1212 

17.6 

82 

90.5 

.120 

17.2 

.75 

6.51 

.388 

1.01 

26.2 

80 

90.5 

.120 

65.2 

1.11 

6.35 

.558 

1.07 

1213 

29 

83.5 

82 

.120 

79 

1.23 

6.64 

.645 

1.10 

1214 

70 

33.3 

41.5 

.078 

— 

2.98 

6.25 

.624 

1235 

16.6 

85 

85.5 

.120 

. 

.71 

6.75 

.378 

17.6 

68 

86.0 

.120 

.75 

5.40 

.321 

.... 

11.7 

66 

85.5 

.120 

51.0 

.79 

5.25 

.425 

1.04 

20.0 

67 

84.0 

.120 

41.8 

.85 

5.32 

.356 

1.04 

21.1 

65 

84.0 

.120 

48.8 

.90 

5.16 

.368 

1.04 

22.9 

63 

83.5 

.120 

46.8 

.97 

5.00 

.384 

1.05 

24.4 

55.5 

84.0 

.120 

51.5 

1.04 

4.41 

.363 

1.03 

25.4 

61 

83.0 

.120 

46.0 

1.08 

4.85 

.414 

1.05 

29.7 

52.5 

83.0 

.120 

31.2 

1.26 

4.17 

.415 

1.04 

31.6 

61.5 

88.0 

.120 

1.34 

4.88 

.517 

34.9 

59.0 

85.5 

.120 

...... 

1.48 

4.68 

.549 

..... 

36.6 

55.5 

84.0 

.120 

. 

1.56 

4.41 

.543 

38.6 

48.5 

82.0 

.120 

. 

1.64 

3.86 

.500 

1237 

26.6 

60.0 

77.5 

.119 

1.13 

4.85 

.*26 

..... 

27.3 

61.0 

78.5 

.119 

1.16 

4.93 

.444 

..... 

21.0 

67.5 

78.5 

.119 

1.19 

5.45 

.505 

..... 

29.6 

57 

79.0 

.119 

_ 

1.26 

4.60 

.452 

31.1 

50.0 

77.5 

.117 

48.6 

1.32 

4.17 

.415 

1.06 

32.3 

53.5 

73.5 

.115 

23.2 

1.37 

4.64 

.461 

1.03 

34.0 

53.2 

78.5 

.119 

-  9.7 

1.44 

4.30 

.481 

.99 

36.6 

48.9 

78.5 

.119 

37.2 

1.56 

3.94 

.460 

1.06 

39.5 

46.2 

80.0 

.120 

32.9 

1.68 

3.66 

.466 

1.06 

43.6 

46.3 

76.0 

117 

... 

28.8 

1.85 

3.87 

.537 

1.05 

FIG.  A-l,  Burning  Rate  for 
Propellant  A-35. 


FIG.  A-2.  Burning  Rate  for 
Propellant  A-146. 


FIG.  A-5.  Burning  Rate  for 
Propellant  A-181. 


FIG.  A-6,  Burning  Rate  for 
Propellant  A-155. 


FIG.  A-7.  Burring  Rate  for 
Propellant  A-156. 


FIG.  A-8.  Burning  Rate 
Propellant  A-157. 


FIG.  A-9.  Burning  Rate  for 
Propellant  A-160. 
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Appendix  B 

INSTABILITY  TESTS  ON  CANDIDATE 
BIG  BOOSTER  PROPELLANTS 


Two  propellants  that  had  been  designed  for  snd  were  candidates  for 
use  in  the  large  solid  booster  program  were  tested  in  both  the  L*-burner 
and  the  low-frequency  T-burner.  They  were  not  tested  in  the  convention¬ 
al  T-burner  because  of  the  high  metal  loading  in  the  propellsnts  which 
would  csuse  high-frequency  oscillations  to  be  completely  attenuated. 

Testing  in  the  L*  burner  wss  more  extensive  than  that  of  the  "big 
pipe"  due  mainly  to  the  ease  of  testing  snd  the  smaller  amount  of  pro¬ 
pellant  used  in  the  L*-burner.  Instability  was  observed  for  both  pro¬ 
pellsnts  with  the  bounds  that  are  normally  observed  in  this  burner.  The 
propellsnts  both  behaved  in  what  would  be  termed  a  "normal"  fashion  and 
on  this  basis  3re  considered  as  being  typical,  i.e.,  nothing  occurred 
that  would  indicate  that  they  would  be  more  prone  to  instability  thsn 
any  other  propellsnt.  The  dats  that  were  obtained  sre  included  in 
Tables  B-l  snd  B-2. 


TABLE  B-l.  Data  Obtained  With  ANB-3105 


Run 

no« 

L*.  in. 

Frequency, 

cps 

Pressure, 

psla 

Growth 

constant, 

sac-1 

1008 

11.4 

59 

80 

<  0 

1009 

7.8 

58 

66 

>  0 

8.7 

53 

70 

'  0 

1010 

13.6 

44 

48 

>  0 

15.2 

46 

50 

>  0 

16.7 

44 

55 

<  0 

1011 

25.1 

39 

52 

'  0 

1012 

25.1 

37 

50 

*  0 

26.2 

37 

57 

'  0 

27,9 

39 

67 

'  0 

1013 

17.8 

43 

48.5 

'  0 

22.3 

45 

56.5 

'  0 

1014 

8.86 

52 

71.0 

19.8 

9.62 

52 

70.5 

'  0 

10.4 

51.6 

70.0 

-12.7 

16.56 

53.0 

83.5 

'  0 
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TABLE  B-2.  Data  Obtained  With  ANB-3254  Propellant 


Run 

no. 

L*,  in. 

Frequency, 

cps 

Pressure , 
psia 

Growth 

constant 

sec"l 

1039 

22.4 

50 

60.5 

70 

21.8 

54 

78.5 

0 

1061 

.... 

87 

.  t  a  . 

1062 

*  •  .  • 

.... 

98 

a  1  a  a 

1070 

35 

35.6 

70.5 

24. b 

.... 

a  •  .  . 

100.5 

a  •  •  a 

1072 

38.5 

30 

50.5 

-14.7 

1073 

47.5 

33 

67 

-13.5 

1081 

27.8 

40 

67.5 

a  »  a  • 

26.4 

37.5 

70.5 

-16.3 

30.7 

31 

62.5 

29.5 

1082 

'35 

41  . 

67 

'  0 

1083 

'69 

30 

57.5 

-10.9 

1084 

7.6 

65 

92.5 

-25.7 

1085 

'12 

55 

68.5 

-21.8 

1086 

-  7.2 

75 

105 

33.0 

1087 

7.6 

75 

95 

-29.0 

1088 

20.9 

50 

73.5 

32.4 

1089 

*27 

43 

83.5 

-18.2 

LOW  FREQUENCY  ACOUSTIC  TESTS 

Seven  tests  were  conducted  in  the  5.5-inch  T-burner  rsnglng  from 
•  15  to  250  cps.  The  pressure  control  system  wss  progrsmmed  to  provide 

for  a  progressively  lncressins  pressure  during  burning  of  the  propellent. 
The  high  burning  rate  of  these  two  propellents  resulted  in  s  tendency 
for  the  relstively  greet  volume  of  combustion  products  to  override  the 
pressure  control  system  snd  to  create  a  burner  pressure-time  history 
that  deviated  considerably  from  that  which  wss  desired.  No  oscillatory 
behavior  was  observed  which  could  be  assigned  to  combustion-driven  pro¬ 
cesses,  except  in  one  test  with  ANB-3254.  Some  mild  oscillations  were 
nhserved  in  other  test?  which  appeared  to  be  caused  by  vibration  of 
pressure  transducers  or  Induced  by  the  pressure  control  apparatus.  Test 
information  is  detailed  in  Table  B-3. 

CONCLUSIONS 

As  a  result  of  testing  these  two  propellants,  it  was  concluded  that 
.  the  possibility  of  either  of  these  propellants  experiencing  pressure 

coupled  acoustic  Instability  or  nonscoustic  instability  in  the  260-inch 
motor  configuration  would  be  rather  remote.  The  possibility  of  velocity 
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coupled  Instability  vaa  not  tested  for  and  could  not  be  ruled  out.  How¬ 
ever,  although  the  poaalblllty  of  thia  type  of  instability  occurring 
exists,  it  vaa  considered  unlikely  In  the  260-inch  configuration. 


TABLE  B-3.  5.5-lnch  T-Burner  Tests 


Propel¬ 

lant 

Test 

no. 

Burner 

length, 

ft 

Pressure 

range, 

psig 

Preaaure 
range 
of  Osc. 
psig 

Freq. 

Remark a 

ANB 

3105 

4820 

12 

20-63 

... 

-*■ 

Low  amplitude  signal 
^  250  cpa  at  end  of 
burn 

4776  i 

24 

24-70 

•  .  . 

Low  amplitude  300  cpa 
signal  above  25  palg 

4777 

l 

24 

0-35 

•  .  . 

-► 

(May  be  transducer 
sctlng  up) 

4821 

6 

0-72 

55-72 

250 

Small  amount  of  high¬ 
er  frequency  present 

ANB 

4819 

12 

40-90 

»  *  « 

Stable 

Poor  preasure  control 

3254 

4778 

24 

30-63 

Stable 

Poor  preasure  control 

4818 

48 

0-55 

Stable 

2nd  mode  In  last  por¬ 
tion  of  teat-probably 
due  to  air  controller 
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